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Abstract 


This survey summarizes briefly results obtained recently in the Casimir en¬ 
ergy studies devoted to the following subjects: i) account of the material 
characteristics of the media in calculations of the vacuum energy (for exam¬ 
ple, Casimir energy of a dilute dielectric ball); ii) application of the spectral 
geometry methods for investigating the vacuum energy of quantized fields 
with the goal to gain some insight, specihcally, in the geometrical origin of 
the divergences that enter the vacuum energy and to develop the relevant 
renormalization procedure; hi) a universal method for calculating the high 
temperature dependence of the Casimir energy in terms of heat kernel coeffi¬ 
cients. A special attention is payed to the mathematical tools applied in this 
field, namely, to the spectral zeta function method and heat kernel technique. 
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I. INTRODUCTION 


In 1948 Casimir [1] proceeding from the general principles of quantum theory of electro¬ 
magnetic held has shown that two uncharged perfectly conducting plates, placed in vacuum, 
should attract each other with the force per unit area 


TT^ch 

240 ^’ 


( 1 . 1 ) 


where a is the distance between the plates. Later this prediction was verihed experimen¬ 
tally [2], For dielectric bodies the Casimir force was also measured with reasonable agreement 
to the theory [3]. 

The Casimir force is very weak, however it increases rapidly as the separation a decreases 
and it becomes measurable when a ~ 1 pm or less. For plates 1 cm in area with a = 0.5 pm, 
the Casimir force is about 0.2 dyn. Firstly the Casimir force has been measured, due to 
technical reasons, between a conducting plane and a sphere with progressively higher pre¬ 
cision using torsion balances [4], atomic force microscopes [5], and capacitance bridges [6]. 
Only recently this force was measured between parallel metallic surfaces in the range 0.5 - 
3.0 pm at the 15% precision level using cantilevers [7]. The Casimir force may be relevant 
in nanotechnology [6,8,9]. At a separation of 10 nm, this force is ~ 1 atm. Now the temper¬ 
ature dependence of the Casimir forces becomes a feasible task for experiment [10,11]. The 
juxtaposition of the experimental results and theoretical calculations with detailed analy¬ 
sis of errors and precision achieved in this held can be found in ref. [11]. An interesting 
proposal was made recently to measure variations of the Casimir energy by making use a 
superconducting cavity [12]. 

The Casimir effect is, in fact, the hrst macroscopic effect predicted in the framework of 
the relativistic quantum held theory, namely, in quantum electrodynamics. The macroscopic 
feature of the Casimir ehect implies that it is observed in the experiments with macroscopic 
bodies, instead of the systems of atomic size. In this respect, it is analogous to supercon¬ 
ductivity and superhuidity. 

This ehect was always an interesting subject of both theoretical and experimental studies. 
Now the Casimir ehect is treated in a more general way, namely, as the inhuence on the 
physical characteristics of a quantum held system of the external conditions (boundary 
conditions, background helds and so on) that lead to the restriction of the conhguration 
space of the system under study. 

Calculations of the Casimir forces for conhgurations more complicated than two parallel 
plates proved to be difficult, and the progress in this held is rather slow. Till now one has 
no intuition even as to whether the force should be attractive or repulsive for any given 
boundary geometry, not to mention the estimation of its value. 

In the case of perfectly conducting surfaces placed in vacuum the Casimir energy is 
calculated exactly, safe for parallel plates [1], for spheres [13-18] and for circular inhnite 
cylinders [19-21]. 

Dielectric properties of the media separated by plane boundaries did not add new math¬ 
ematical difficulties [22-24] as compared with the Casimir pioneer paper [1]. However the 
hrst result on the calculation of the Casimir energy for the non-hat boundaries was obtained 
only in 1968. By computer calculations, lasted 3 years, Boyer found the Casimir energy of 
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a perfectly conducting spherical shell [13]. Account of dielectric and magnetic properties of 
the media in calculations of the vacuum energy for nonflat interface leads to new difficulties 
or, more precisely, to a new structure of divergencies. 

The calculation of the Casimir energy in a special case, when both the material media 
have the same velocity of light, turns out to be, from the mathematical standpoint, exactly 
the same as for perfectly conducting shells placed in vacuum and having the shape of the 
interface between these media [19,25-36]. 

The hrst attempt to calculate the Casimir energy of a dielectric compact ball has been 
undertaken by Milton in 1980 [37]. And only just recently the hnal result was obtained for 
a dilute dielectric ball at zero [38-45] and hnite [46,47] temperature. 

Situation with a dielectric compact cylinder proves to be much more complicated. The 
main reason of this is that here there is no decoupling between transverse-electric and 
transverse-magnetic modes [19]. It was shown that the summation of van der Waals in¬ 
teractions inside a purely dielectric cylinder in the dilute-dielectric approximation lead to a 
surprising null result [19]. First it was derived independently by Romeo and Milonni (un¬ 
published calculations) and by Barton in ref. [48]. Bordag and Pirozhenko [49] have found 
the hrst few heat kernel coefficients for this conhguration. Their results imply, in particular, 
that zeta function technique should supply a hnite value for vacuum energy of a dielectric 
cylinder in the dilute approximation. And only very recently Cavero-Pelaez and Milton [50] 
have shown rigorously (in quantum held theory approach) that the Casimir energy of an 
inhnite circular dielectric cylinder vanishes through the second order in the deviation of 
the permittivity from its vacuum value. This result is of a hve-year-long ehort! The same 
vanishing value for this energy was derived anew in paper [51] by making use of the mode 
summation method. 

For the force between a plane and a sphere the exact result is not known, and only an 
estimation valid when both are close enough is available. This estimation is based on the 
so called proximity theorem [52,53]. Similar estimations have been recently obtained for 
the force between two spheres [54]. The problems of two concentric spheres [45,55] and two 
concentric cylinders [56] have also been considered recently. For both cases expressions for 
the Casimir energy have been obtained using the Abel-Plana formalism. However a detailed 
numerical calculation and analysis of the results are still missing. The authors of ref. [57] 
studied the Casimir interactions, due to the massless scalar field fluctuations, of two surfaces 
which are close to each other. Two close co-axial cylinders, co-centric spheres co-axial cones 
and two tori were considered. The first correction to the parallel plates result was calculated. 

For dielectric-diamagnetic media with a continuous velocity of light at the interface the 
exact results are known for a sphere [25,27,29-36] and cylinder [19,26,35]. 

For arbitrary dielectric media the exact results are obtained for parallel plates (Lifshitz 
formula [22-24]), for a dilute dielectric ball [38-43,46] and for a dilute dielectric cylin¬ 
der [50,51]. A single surface immersed in an inhomogeneous medium also experiences the 
Casimir force [58]. 

Recently there appear comprehensive reviews [59,60] and excellent book [61] which are 
devoted to the Casimir effect. Previous survey publications in this held [62-67] remain 
valuable. Many references concerning vacuum energy can be found in [68]. Of special 
interest are the Proceedings of regularly organized, every 3 years since 1989, Workshops on 
Quantum Field Theory Under the Inhuence of External Conditions [69-72]. 
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In this situation a new survey in this field should be naturally concerned with particular 
problems to investigation of which a certain contribution has been done by present authors. 
First of all, we are going to give a lucid and at the same time comprehensive introduction 
into the mathematical methods that are used in current Casimir calculations. We hope that 
the review will provide a sufficient material in order to apply these methods to coinciding 
problems in diverse branches of theoretical physics. 

The layout of the review is as follows. In sect. II we consider the basic terms used in 
the Casimir studies, specifically, the notion of zero point energy. To our mind employment 
of this term is not only convenient but also physically justified in this field. In sect. Ill the 
mathematical methods applied in Casimir calculations are considered, namely, the spectral 
zeta functions and heat kernel technique. In sect. IV the spectral zeta functions for sim¬ 
plest boundary configurations are calculated and the relevant Casimir energies are found. 
Section V is devoted to calculation of the vacuum energy of electromagnetic field connected 
with a dielectric ball. This problem proved to be nontrivial one. We are trying to reveal 
the principal drawbacks encountered here and summarize the experience obtained in these 
calculations. The physical origin of the divergencies in vacuum energy of quantized field 
and their relation to the boundary properties in considered in sect. VI. In sect. VII a uni¬ 
versal method of constructing the high temperature asymptotics in the Casimir calculations 
is discussed. It is substantially based on the heat kernel technique, more precisely, it uses 
the coefficients of the asymptotic expansion of the heat kernel. In sect. VIII an attempt is 
undertaken to summarize the lessons gained in the course of quiet long studies in this field 
and to formulate the present status of our physical understanding of the Casimir effect. In 
the Appendices A, B and C the mathematical details of the calculations involved are placed. 


II. ZERO POINT ENERGY IN QUANTUM MECHANICS AND IN QUANTUM 

FIELD THEORY 

The quantum theory forbids the particle to posses simultaneously a definite value of its 
coordinate and momentum. This results, in particular, in the following. The particle placed 
in a smooth potential with a local minimum cannot acquire the value of energy which 
corresponds to the potential minimum. In fact, smooth potential can be approximated 
around its minimum by the potential of the harmonic oscillator with the known spectrum 

En = lTuj(^+^, n = 0,l,2, - (2.1) 

Thus the particle cannot posses the energy less than Eq = fvjj/2, which has obtained the 
name of zero point energy. It is interesting that the notion of the zero point energy has first 
appeared before the development of quantum mechanics, in Planck paper of 1912, where he 
has suggested the “second theory” of the black body spectrum [73]. By making use of zero 
point energy, Einstein and Stern derived the Planck distribution following practically the 
classical consideration [74]. 

The prediction of zero point energy, obtained in the framework of the quantum mechanics 
dealing with the systems of a finite number of degrees of freedom, is directly verified in 
experiment. For example, by making use of the Roentgen rays or neutrons scattering by a 
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crystal lattice or by liquid helium one can convince that, when temperature tends to zero, 
the atoms in these systems indeed occupy the state with the energy Eq and they cannot 
possess lower energy. In actual fact these experiments provide the Debye-Waller factor [75], 
calculation of which in the framework of the quantum mechanics takes into account the zero 
point oscillations of the atoms. 

A direct justihcation of the reality of zero point oscillations (or zero point energy), in the 
case of a quantum mechanical systems, is afforded by the observation of vibration spectra 
of diatomic molecules. It is interesting that it has been also done prior to the development 
of quantum mechanics [76]. 

The zero point energy is also manifested in the temperature dependence of the specihc 
heat of rigid bodies. For example, in the model suggested by Einstein [77] the crystal lattice 
is considered as the system of 3iV independent oscillators, N being the total number of lattice 
sites. All the oscillators have the same frequency (a single phenomenological parameter of 
the model). In this model the area between the curve c(T) describing the temperature 
dependence of the specific heat and the straight line y = c(cx)) is proportional to the zero 
point energy of all the oscillators. In a more realistic model of specihc heat of crystal bodies 
proposed by Debye [78] the zero point energy is also taken into account without problem. 

Drastically different situation concerning the zero point energy arises when one considers 
the system of infinite number of oscillators. It is this case that is encountered in quantum 
held theory. The Hamiltonian of the free (noninteracting) quantum held is an inhnite snm 
of the Hamiltonians of the harmonic oscillators. When the held is considered in unbounded 
space this sum is substituted by the integral 

H = - J (i^kci;(k)[a^(k)a(k) -|- a(k)a^(k)], (2.2) 

where the frequency Cc;(k) is equal to the energy of the particle with momentum k. For 
simplicity we are considering uncharged spinless particles (bosons) 

[a(k), a^(k')] = (5^^^(k — k'). (2.3) 

The eqnidistant energy spectrum of the harmonic oscillator (2.1) enables one to interpret the 
qnantnm held dynamics in terms of particles, namely, the transition of one of the oscillators 
in the Hamiltonian (2.2) to the next higher (lower) level is treated as creation (annihilation) 
of the particle with the energy a;(k). In the framework of this interpretation, the zero point 
energy of qnantnm held Eq is treated as the vacuum energy, because in this case there are no 
excited oscillators (all them occnpy the energy level Eq). Hence there are no real particles. 
For any physically acceptable dispersion law Cc;(k) this energy is, obviously, unbounded 



Naturally, question arises here how to treat this energy. 

Divergences are typical for any realistic quantum held theory because the Hamiltonian 
of qnantnm held (see, for example, the free held Hamiltonian (2.2)) involves the oscillators 
with arbitrary large freqnencies ci;(k) when k ^ oo. Physically it is clear that such quanta 
cannot be relevant in any concrete calculation. However, qnantnm held theory does not 
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provide us with the mechanism that could suppress the contribution of such quanta, and the 
problem of divergencies in relativistic quantum field theory requires suitable renormalization 
procedures. 

If the energy in the system under consideration is conserved (for example, the scattering 
of elementary particles possessing the same total energy in the initial and hnal states is 
considered) then the inhnite vacuum energy Eq can be ignored by shifting the origin of the 
energy scale by Eq. Formally this subtraction is accomplished by transition to the normal 
ordering of the operators in the held Hamiltonian (2.2) 


H = 


1 

2 


(i^kci;(k)[: a^(k)a(k) : + : a(k)a^(k) :] 


J d^k ct;(k) a^(k)a(k). 


(2.5) 


However when the energy of the quantum held is not conserved then such a procedure 
of removing Eq cannot be justihed. This issue is of a special importance in the gravitation 
theory because the energy of matter helds is the source of the gravitation held. 

In quantum held theory there is only one receipt of a consistent, from the mathematical 
standpoint, treatment of the divergencies. It is the renormalization procedure, i.e., transi¬ 
tion from the initial parameters of the theory, which turn out to be, as a rule, inhnite and 
therefore unobservable, to hnite physical parameters [79,80]. On principle, the renormaliza¬ 
tion procedure in quantum held theory should be done always independently of the presence 
of divergencies. If the theory does have the inhnities then they are removed in the course of 
renormalization “by the way”. Certainly, it is true only for renormalizable held theories. 

A consistent renormalization procedure is developed only for the scattering processes in 
unbounded Minkowski space (for example, the Bogoliubov R operation). If the held occupies 
a bounded region in space, the relevant renormalization scheme is absent. For example, the 
renormalizability of quantum electrodynamics considered in a compact space region, by 
imposing the corresponding boundary conditions, is not proved explicitly. However there 
is no reason to suspect this renormalizability. Indeed, the boundary conditions, conhning 
the helds to a compact region, arise due to transition from a microscopic description of the 
problem at hand (in terms of charges and electromagnetic held by making use of standard 
renormalizable quantum electrodynamics in unbounded space) to a macroscopic description 
by imposing boundary conditions on electromagnetic held and introducing phenomenological 
characteristics (permittivity and permeability) specifying the materials the boundary made 
up. 

With lacking general renormalization algorithm for nontrivial boundary conditions, one 
has to develop this procedure for each concrete conhguration a new. In fact, it is the main 
drawback in the Casimir calculations for boundaries of diverse geometries. Undoubtedly, 
comparison of the theoretical predictions with the experimental observations should play 
here the crucial part. By now only the simplest conhguration, plane plates, is feasible for 
experimental study (see Introduction). In the next sect, we consider diherent mathematical 
methods for removing the divergencies in the Casimir calculations. 

It is interesting to note that in the case of fermions the vacuum energy has opposite sign 
in comparison with bosons. Really, the Hamiltonian function for noninteracting fermionic 
oscillators is given by [81] 

Hf = ^ [ (i^ka;(k)[&^(k)&(k) — 6(k)6'^(k)]. (2.6) 


6 



In quantum theory the amplitudes 6(k) and b^(k) obey the anticommutation relation 

|6(k),M(k)]+ = iW(k-k'), (2,7) 


which enables one to cast eq. (2.6) to the form 

Hf = j d^ko^(k) b\k)b{k) + Eof, (2.8) 

where Eqf is the fermionic zero point energy 

EoF = ~j d^kuik). (2.9) 

As known, the opposite signs in eqs. (2.4) and (2.9) result in substantial conciliation of 
divergencies in held theories symmetric with respect to permutation of bosons and fermions, 
i.e., in supersymmetric theories. 

The Casimir effect can be considered without introducing the notion of the vacuum 
helds and its zero-point energy, for example, in the framework of the Schwinger’s source 
theory [82-84]. In this connection one can meet the opinion [65], according to which there 
is no sense to insist that “... observable phenomena like the Casimir effect strongly suggest 
that the vacuum electromagnetic held and its zero-point energy are real physical entities 
and not mere artihces of the quantum formalism”. And it is largely “matter of taste” to use 
for the interpretation of this ehect vacuum or source helds. 

A similar assertion has been done in a recent paper [85] by R. L. Jahe: “... no known 
phenomenon, including the Casimir ehect, demonstrates that zero point energies are “real”. 
The author stressed once more the well known fact that the Casimir force originates in the 
forces between charged particles in the metal plates or, more precisely, between huctuating 
dipoles inside the diherent plates [22,61,86]. In view of this, it is obvious that in the micro¬ 
scopic theory, describing the atomic structure of the plates, the Casimir force vanishes as 
a, the hne structure constant, goes to zero. The independent of a result (1.1) obtained in 
macroscopic approach corresponds to the limit a —cx) in the microscopic consideration. It 
is clear that in microscopic description the notion of zero point energy does not appear. 

However we think that denying the reality of the zero point energy of physical fields 
mentioned above is not sufficiently justihed. Here it is worthy to recall that the usual classical 
mechanics may be reformulated without using the notion of the force, as has been proposed 
by Hertz [87] (see the evaluation of this approach due to Sommerfeld [88]). However, now 
nobody suspects the reality of forces, as the physical reasons that cause the alteration of 
the velocity of a test particle. Therefore we believe that it is well-grounded to treat the 
experimentally observed Casimir forces as the direct manifestation of the zero point energy 
of the relevant vacuum helds [61] because the use of this term enables one to accomplish the 
theoretical analysis of the Casimir effect in the most simple and clear way. 

Closing this section it is interesting to note that the notion of the zero point energy was 
also considered and used in the framework of the classical theory of electromagnetic held 
(see the papers [89-91] and references therein). 
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III. RENORMALIZATION OF THE VACUUM ENERGY 


A. Physical background for removing the divergencies 

In the general case one cannot assign the physical meaning to the vacuum energy of 
quantum held (2.4) unlike to the zero point energy of quantum mechanical systems with 
hnite degrees of freedom. However if the conhguration space M is bounded (the Casimir 
effect), then one can get a hnite value for the vacuum energy subtracting appropriate counter 
terms, i.e., by making the renormalization. 

Certainly, in the experiment one can observe only the Casimir forces that are derived 
from the vacuum energy by diherentiation with respect to relevant distance. We shall always 
bare in mind this when considering the vacuum energy problem. 

When the held occupies a bounded region M, its eigenfrequencies are discrete and their 
values are determined by the geometry of the boundary dAi. As a result, the vacuum energy 
also depends on the form of the boundary 

Bo(dM) = ^^u;n. (3.1) 

n 

Following the basic line of the quantum held theory formalism, one can expect that an 
observable value of the vacuum energy is obtained by subtracting from eq. (3.1) the same 
expression calculated for the boundary of a special form. In the simplest case, when elec¬ 
tromagnetic held is conhned between two parallel perfectly conducting plates, it is sufficient 
to subtract from eq. (3.1) the contribution of unbounded Minkowski space i?o(0) 

Eo = EoidM) - Eo{0). (3.2) 

Preliminary both expressions for Eq should be regularized, for example, by introducing a 
smooth function, that suppresses the contributions due to large momenta. In this way the 
Casimir energy for parallel plates is calculated, for example, in the text book on quantum 
held theory [92]. 

The prescription (3.2) for obtaining the physical value of the vacuum energy may be in¬ 
terpreted in the following way. In the case of inhnite number of degrees of freedom (quantum 
held theory) the observable quantity is not the zero point energy itself, but only its excess, 
caused by boundaries or inhomogeneities of the space, compared with the zero point energy 
of quantum held occupying unbounded homogeneous space. It is the main point that dihers 
treating the zero point energy in quantum mechanics and in quantum held theory. 

However, the subtraction of the Minkowski space contribution turns out to be insufficient 
in order to obtain hnite vacuum energy in the case of more complicated boundaries, for 
example, for perfectly conducting sphere. Unfortunately, there are no general rules that 
could enable one to construct the subtraction procedure for the boundary of arbitrary form, 
i.e., , to hnd the corresponding counter terms. In view of this, in calculation of vacuum 
energy the analytical regularization, or zeta regularization, proved to be very useful because 
it simultaneously accomplishes the renormalization without introducing the relevant counter 
terms. Here we are dealing with renormalization of the held energy. This problem is central 
one in quantum held theory on nonhat background [93]. 



There is a few surveys and books concerned with the spectral zeta function technique and 
its application to the Casimir calculations [60,94-102], Therefore we present here only the 
basic formulae and stress the points that are not sufficiently elucidated in the mathematical 
literature. 


B. Spectral zeta functions 

Let the held (p(t,x) is governed by the following equation 

+ = 0 ’ ( 3 - 3 ) 

where L is an elliptic differential operator of the second order acting only on the spatial 
variables x. For example, in the case of a free scalar massless held L = —A, where A is the 
Laplace operator. 

Upon separating the time dependence (p(t,x) = the held equation (3.3) and 

the relevant boundary conditions with respect the variable x generate the spectral problem 

L(fn{x) = Xn<fn{x) or L\n>= Xn\n>, Xn = uil/c^. (3.4) 

It is convenient to use here the Dirac bracket notation. We assume that the spectral problem 
is well posed, i.e., , all the equations written below have sense. 

The completeness relation of the vector set {|n >} can be represented (formally) in the 
following way 




n >< n 


n 


(3.5) 


where / is a unity operator acting in the linear space of the vectors \n >. In view of this we 
have for the inverse operator L~^ 


L-^ 


E 


n >< n 


X 


n 


(3.6) 


It can be easily checked with allowance for eq. (3.4). In the same way we have for the sth 
power of the inverse operator L~^ 




E 


n >< n 


A 


s 

n 


(3.7) 


The local spectral zeta function (l{s]x) of the operator L is a diagonal element of the 
operator 

A / \ ^ x\n X n\x > \ x / x /o o\ 

Cl{s-, a:) = -= zZ <^3.8) 

n n 

The global spectral zeta function [103,104] is dehned by 
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(3.9) 


a(s)=Tl'L-* = 5^Ar. 

n 

It is obtained by integration of Cl('S; x) over the whole space 

Cl{s) = j CL{s-,x)dx. (3.10) 

The dehnition of the spectral zeta function (3.9) is a direct extension of the Riemann 
zeta function [105] 

oo 

Cr(^) = X]—) Res>l (3.11) 

n=l 


to the spectrum of the operator L. 

In the case of the d-dimensional space the index n in eq. (3.9) stands for the set of d 
indices n = {ni,n 2 , ■ ■ ■ ,nd} and summation should be carried out with respect to each of 
them in respective ranges. In the majority of the problems, having physical application, 

Xn ^ nl + nl + ... + ri^, n* —> cx), 1 <i < d. (3.12) 

In this case the series (3.9), dehning the global zeta function, converges in the domain 
Re s > d/2, where it represents an analytical function Cl{s) of a complex variable s. This 
function can be analytically continued into the left half-plane Re s < d/2 except for simple 
separate poles, for example, at the point s = d/2. This fact is a mathematical basis for 
employment of the spectral zeta functions with the aim of obtaining hnite expressions when 
treating the divergencies. 

In the framework of the zeta function technique the vacuum energy Eq of the quantum 
held (p{t,x) is dehned by [99] 



In view of eqs. (3.10) and (3.13), the quantity {1/2)(l{—1/2]x) can be interpreted as the 
vacuum energy density. Besides this, it is easy to show that the vacuum expectation value 
of the canonical energy-momentum tensor Tqq{x) for the held (p{t,x) can be represented in 
the form 

< 0|Too(a;)|0 >= ]^'^Xll^^l{x)^n{x) = • (3.14) 

In order to construct the global zeta function, the spectrum of the problem under consid¬ 
eration should be known. Certainly, it is feasible only for boundaries with high symmetry, 
for example, with spherical or cylindrical symmetry. In the Casimir calculations with d = 3 
and n = {ni,n 2 ,n 3 }, the conhgurations are usually considered when, for hxed values of rii 
and n 2 , the eigenfrequencies are the roots of the equation 

fn,n,{oo) = 0 (3.15) 
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with a known function fnin 2 iz)- The sum of these roots can be found by making use of 
a contour integration in complex frequency plane [105]. It leads to the following compact 
representation for the spectral zeta function 

= (3-16) 

ni,n2 ^ 

where the contour C encloses, counterclockwise, all the positive roots of eq. (3.15). This 
representation proves to be convenient for analytical continuation of the right hand side 
of eq. (3.16) to the left half plane, specihcally, at the point s = —1/2. However, it may 
happens that the zeta function (/(s) has a pole at this point. It will imply that the zeta 
regularization does not afford a hnite value for the vacuum in the problem at hand. There 
is a rigorous mathematical criterion that enables one to predict when this happens. In order 
to formulate it we need another spectral function, heat kernel. 

Before to go over to consideration of this function it is worthy to make a general comment 
concerning the spectral functions. In the theory of differential operators these functions 
play in fact the same part as characteristic polynomials do for matrix operators of a hnite 
dimension, namely, they store all the information about the spectrum of a given operator in 
an “encoded” form. When the spectrum consists of a hnite number of eigenvalues A„, n = 
1, 2,..., iV , it is natural to construct on this basis the polynomial of iVth degree the roots 
of which give this spectrum; PnW = 11^=1 (-^ ~ ^n), A G C. The matrix equation 
P]\f{Lj^) = 0 is used, as known, for constructing the functions of the matrix L^r. In the case 
of an inhnite spectrum a concrete form of the relevant spectral function is not so dehnite. 
Ultimately, all is determined by the problem under study. For example, when dealing with 
the spectral sums, it is convenient to use the zeta function (3.9). 

The procedure of constructing the local spectral zeta function (3.8) proves to be much 
more complicated as compared with that for the global zeta function. Besides the spectrum 
of the operator L we have to use the respective natural modes. As far as we know the local 
zeta function has been derived only for the scalar Laplace operator dehned on a cone or in 
wedge. This function was used in studies of black hole physics [106,107], cosmic strings [108] 
and in calculation of the Casimir effect for wedge [109]. 

The relation of the zeta function technique to more transparent regularizations is dis¬ 
cussed, for example, in refs. [110,111]. 


C. Heat kernel technique 

In Casimir calculations also useful is another spectral function of the operator L intro¬ 
duced in eq. (3.3), namely, the heat kernel of this operator 

K{t) = Tr (e-"^) = ^ 6"^"^ (3.17) 

n 

where r is an auxiliary variable ranging from 0 to -|-cx). Such a name of this function is due 
to the following. By making use of the unity operator (3.5) one can write 
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(3.18) 


= 


E- 


,-tA„ 


\n >< n\. 


The matrix element of this operator 
K{x,y]T) = < x\e-^^\y > 


< x\n >< n\y >=^e 


is the Green function of the heat conduction equation with the operator L 

d \ 

Lx + K{x,y]T) = 0 , 


(3.19) 


(3.20) 


y; r) =5{x,y), r ^ +0. 

For the functions K{t) and K{x,y; r) the relation analogous to (3.10) holds 

A'(r) = f dx K{x, x; t). 


(3.21) 


(3.22) 


The integrated heat kernel (3.17) and the spectral zeta function (3.9) are connected by 
the Mellin transform. Indeed, from the dehnition of the gamma function [112] it follows that 


Af 


r(s) 


dTT‘-‘e 


S —1 „ — 


Re s > 0. 


(3.23) 


Upon summation in both sides of this equation we obtain 


cw = 


r(s) 


dxT^ K{t), Res> 


d 


(3.24) 


To the left half-plane Res < (i/2 eq. (3.24) should be continued analytically. 

Thus, from the mathematical point of view, both the spectral functions C(>s) and K{t) 
are equivalent. However in practical use it is difficult to succeed in hnding the spectral zeta 
function or the heat kernel in an exact form, and one has to deal with the approximate 
expressions. Therefore different initial formulae either (3.9) or (3.17) and (3.20) - (3.22) 
prove to be useful. 

In physical applications are important the coefficients in the asymptotic expansion of the 
heat kernel, when r —-|-0 


K{t) = = (47rr)-'^/2 ^ + ES. (3.25) 

n n=0,l,2,... 

In this expansion d is the dimension of the configuration space manifold M. in the problem at 
hand, ES stands for the exponentially small corrections as r —-|-0. The first two coefficients 
in expansion (3.25) are determined by the volume V of the manifold Ai and by the area of 
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the boundary DM.. For example, if L = —A, where A is the Laplacian, acting on a massless 
scalar held, then 


B, = V, By 2 = T^S. (3.26) 

The upper sign in this equation is for the Dirichlet boundary conditions and the lower sign 
is for the Neumann conditions. 

For hat manifolds all the coefficients Bn/ 2 , n>l are due to the boundary contributions, 
namely, they are dehned by the integrals over the boundary with the integrands expressed 
in terms of geometrical invariants of the boundary [60,94,95,97]. The boundary is expected 
to be sufficiently smooth. 

We give here the properties of the heat kernel coefficients which are especially important 
when calculating the vacuum energy. The proof of these statements can be found in the 
literature (see, for example [60,94,95]). The hrst few coefficients Bn /2 yield the ultraviolet 
divergencies of the vacuum energy (3.13). For d = 3 it is the coefficients up to i? 2 - If 
the coefficients B 2 is equal to zero, then the zeta regularization gives a hnite value for the 
vacuum energy according to the formula (3.13). For d = 2 the role of such an indicator 
plays the coefficient B 3/2 (see below). The heat kernel coefficients also specify conformal 
anomalies taking place in a concrete held theory model, the high temperature behavior of 
the thermodynamic functions (see sect. VII) and so on [94]. 

The heat kernel coefficients are determined by the residua of the product F(s)C(<s) at cor¬ 
responding points. Indeed, substituting the expansion (3.25) into eq. (3.24) and integrating 
over r in the vicinity of the origin we obtain 


Bn/2 

(47r)V2 


lim 

^ d—n 

~T 


n — d 


2 


r(s)C(s), n = 0,l,2. 


(3.27) 


In order to use this formula the spectral zeta function (C(s) should be known in the vicinity 
of the following points 


d d — 1 
^ ^ 2 ’ 2 


(3.28) 


For this one has to do analytical continuation of the initial formula (3.9) to the left half-plane 
Re s < d/2. 

If d = 3 and R 2 7 ^ 0 or for a given d the coefficient B{d+i )/2 does not vanish then according 
to eq. (3.27) the zeta function has a pole at the point s = — 1 / 2 . As a result in these cases 
the renormalization by zeta function does not provide a finite value for the vacuum energy 
Eq (see the dehnition (3.13)). 

There is another technique for calculating the heat kernel coefficients that enables one to 
avoid the procedure of analytical continuation [113-116]. It is close to the method presented 
above but sometimes calculations become simpler. 

Let us consider the spectral zeta function depending on a parameter 


C{.s,x^) = ^(A„ + a;2) L 

n 


(3.29) 
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It may be regarded as an extension of the Epstein-Hurwitz zeta function 


CEH(-s,a^) = + 0 ?) " 

n=l 


to the general spectral problem (3.4). 

It turns out that the heat kernel coefficients Bn /2 can be found from the expansion of 
the function ({s,x‘^) in terms of inverse powers of x developed for a certain value of s. It 
is convenient to chose this value to be equal to 1 + (i/2. In fact, from the definition of the 
gamma function (3.23) it follows that 


r(s)(A + a;2)-* = / Res>0. 

Jo 

For s = 1 + d/2 eq. (3.30) gives 


/ CXD 

/ OO 


> — 


On substituting the asymptotic expansion (3.25) in eq. (3.31) we obtain 


d 


d 


r 1 + ^ c in.xd-E 


B. 


n/2 


n=0 


n 


^ (47r)^/2^ 2.'"^ " ' 


(3.30) 


(3.31) 


Bo 5i/2r(3/2) B,T{2) Bo/2r{5/2) 

o A c 


{inyB ' x^ ' x"^ 

, fi2r(3) , 55/2r(7/2) R3r(4) 


X'-^ 


+ 


+ 


X'^ 


x‘ 


x^ 


+ 0{x 


- 9 \ 


(3.32) 


Let F{z) =0 be the frequency equation which determines the spectrum in the problem 
under consideration. We also suppose that the function F{z) allows one to rewrite this 
equation in the form 


JJ(A„ - z‘^) = 0. 


Taking into account the equality 


1 

{Xn + x^)^ 


(_l)m 

r(m) 


d A 
2xdx J 


ln(A„ + x^), 


z = ix, 


we recast the left-hand side of eq. (3.32) to the form 


r 



+ 5 .^ 



2x dx J 


l+d/2 

In F{ix). 


(3.33) 


(3.34) 


(3.35) 
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Obviously formula (3.35) is applicable only to the manifolds of even dimension. Now we 
have to substitute the left-hand side of eq. (3.32) by the expansion of eq. (3.35) for large x^. 
After that the heat kernel coefficients Bn /2 are obtained by comparing the coefficients of 
in both sides of eq. (3.32). 

This technique can be extended to the odd values of d by calculating the zeta function 
({s,x^) at the point s = (1 -|- d)/2. Instead of eqs. (3.32) and (3.35) we have now 


r 




(3.36) 


Again comparing the coefficients of x in both sides of eq. (3.36) we obtain the heat kernel 
coefficients for odd d. 


IV. SPECTRAL ZETA FUNCTIONS AND VACUUM ENERGY FOR SIMPLEST 

BOUNDARIES 

Practically every problem on calculation of the Casimir energy (or force) has been con¬ 
sidered often with employment of more and more effective and elaborated mathematical 
methods. For example, the hrst calculation of the Casimir energy of a perfectly conducting 
spherical shell carried out by T.H. Boyer in 1968 [13] has required computer calculations 
during 3 years [117]. Later this problem was considered in many papers [14-16]. By making 
use of the modern methods [118] it can be solved almost without numerical calculations 
(with a precision of a few percent). It requires only the application of the uniform asymp¬ 
totic expansion for the Bessel functions. The zeta function regularization enables one to 
calculate the vacuum energy in a consistent way without dealing with the manifestly inhnite 
expressions. 


A. Parallel perfectly conducting plates 


As known, for example, from the theory of waveguides and resonators [119] the vectors 
of electric and magnetic helds in the problem at hand are expressed in terms of the electric 
(n') and magnetic (n") Hertz vectors, each having only one nonzero component H^ and 
satisfying, respectively, Dirichlet and Neumann conditions on the internal surface of the 
plates. The functions H^ and obey the equations 


^2 

-h V 

dz^ 


n, = —n„ 


dz‘^ 


+ 


„// UJ ,, 

n, = —n„ 


(4.1) 


where uj is the frequency of electromagnetic oscillations, stands for the two-dimensional 
Laplace operator for the variables (x, y) = x. The separation of variables results in the 
following solution 
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m^z 


n= 1, 2,..., 


nl(x,^) 

n"(x,^) 

^n(k) 


= exp(ikx) sin 
= exp(ikx) cos 


a 

mrz 


= c 




n = 0, 1, 2,..., 


(4.2) 


where a is the distance between the plates. Hence, the states of electromagnetic held with 
the energy tvjJm n > 1, are doubly degenerate, while the state with the energy hojQ = hck 
is nondegenerate. 

With allowance for this the zeta function in the problem under consideration is given by 


C(^) 


Lx Ly 

n2s 


d^k 

(27r)2 




n=l k 


k^ + 



+ (k^ + 



(4.3) 


where and Ly are the dimensions of the plates. 

For a correct dehnition of the integral in this formula in the small k region the photon 
mass /i is introduced (infrared regularization). At the hnal step of calculations one should 
put /i = 0. On integrating in eq. (4.3) and substituting the sum over n by the Riemann zeta 
function one arrives at the result 


C(^) 


LxLy \ CR(2g - 2) 

27rc^^ Va/ s — 1 2s — 1 


(4.4) 


The zeta function (4.4) leads to the well-known value for the Casimir energy 


or for its density 



TT LxLy 


V 


cHtt^ 

720a^’ 


where V 


d L 


y 


(4.5) 


(4.6) 


Differentiation of the vacuum energy (4.5) with respect to the distance a gives the Casimir 
force (1.1). 


B. Sphere 

Let us consider a solid ball of radius a, consisting of a material which is characterized 
by permittivity ei and permeability pi. The ball is assumed to be placed in an inhnite 
medium with permittivity 62 and permeability fi 2 - In the case of spherical symmetry the 
solutions to Maxwell equations are expressed in terms of two scalar Debye potentials ip (see, 
for example, [120]) 

E™ = V X V X (r^™), H™ = x (r^™) (E-modes), 

E™ = tuV X (r^,,.), = V X V X (r^^) (E-modes). (4.7) 
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These potentials obey the Helmholtz equation, have the indicated angular dependence, and 
are regular at the origin 

^2 

(V^ + k‘f)tlJim = 0, i = l,2, {r ^ a); ^|Jlm{r) = (j){r)Yim{kl). (4.8) 

Here we have to make a note concerning the formulation of the spectral problem for 
differential elliptic operators dehned on unbounded manifolds. First of all, the question 
arises what conditions should be imposed on the eigenfunctions at the spatial inhnity. In 
order for the uniqueness theorem to hold for the solutions looked for, the radiation conditions 
at infinity should be introduced. It is a well known fact in classical mathematical physics (see, 
for example, [121]). However, in this case there are no solutions with real frequencies and 
one has to consider the complex frequencies and respective natural modes, i.e., quasi-normal 
modes. Certainly, it is absolutely unclear how to use these modes when quantizing the helds. 
Here one can proceed in two different ways. The quantum held under consideration can be 
placed inside a sphere of a large radius R imposing on held functions reasonable conditions 
on the internal surface of the sphere [18]. All outside the sphere is ignored. The spectrum 
of oscillations in this case is discrete, inhnite and real because we are dealing simply with 
standing waves. On this basis the spectral zeta function can be constructed in a standard 
way and ultimately the radius of an auxiliary sphere R should be allowed to go to inhnity. 

More appealing, from the physical point of view, is an approach that uses at inhnity the 
conditions on held functions taken from the scattering problem instead of the radiation condi¬ 
tions. The frequency spectrum in this case is real, positive but continuous. The drawbacks 
due to unnormalizable natural modes can be overcome by introducing into consideration 
the relevant spectral density. This density is expressed in terms of the phase shifts, or more 
precisely, through the dost function of the corresponding scattering problem. 

Both these approaches lead to the same hnal results. Furthermore, it turns out, that 
one can also deal with the radiation conditions at inhnity and consequently with complex 
frequencies, ignoring the problem of treating the corresponding quasi-normal modes. Point 
is that the equation for real frequencies in the case with additional large sphere reduces, 
when i? —cx), to the frequency equation in the approach which uses the radiation condition 
at inhnity. And further, the dost function mentioned above is determined by the left-hand 
side of the equation for complex frequencies. 

In view of all this, we shall impose at inhnite the radiation conditions on the function 
0(r) in eq. (4.8) and deal with the frequency equation with complex roots keeping in mind 
that it is in fact the relevant dost function. 

At the ball surface the tangential components of electric and magnetic helds are con¬ 
tinuous. As a result, the eigenfrequencies of electromagnetic held for this conhguration are 
determined by the frequency equation for the TE-modes 

Aj^{au) = s'i{kia) ei{k2a) - y/efjll sfkia) e[(A;2a) = 0 , (4.9) 

and the analogous equation for the TM-modes 

A™{auj) = ^/e 2 jHs'i{kia)ei{k 2 a) - ^/eip 2 sfkia) e'i{k 2 a) = 0, (4.10) 

where ki = ^JeffTiOJ, i = 1,2 are the wave numbers inside and outside the ball, respec¬ 
tively [122]. Here sfx) and efx) are the Riccati-Bessel functions 
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Slix) = Ji+i/ 2 {x ), ei{x) = , (4.11) 

and prime stands for the differentiation with respect to their arguments, kia or k 2 a. The 
orbital momentum I in eqs. (4.9) and (4.10) assumes the values 1,2,.... The roots of the 
frequency equations (4.9) and (4.10), including the complex ones, have been investigated by 
Debye in his PhD thesis concerned with the light pressure on a material ball [123]. 

The material presented below is based in main part on the paper [124]. When considering 
the boundaries with spherical and cylindrical symmetries, it is convenient to subtract in the 
dehnition of the spectral zeta function (3.9) the contribution of the same conhguration under 
the condition a ^ oo with frequencies Up 

CW = - i.;*). (4.12) 

P 

As usual when one is dealing with an analytic continuation, it is useful to represent the 
sum (4.12) in terms of the contour integral by making use of eq. (3.16) 


Cball(-S) 


2 / + 1 / dz d Aj^{az)A™{az) 

27 ri p^o Jq {z'^ + dz A™(cx))A™(cx)) ’ 


(4.13) 


where the contour C surrounds, counterclockwise, the roots of the frequency equations (4.9) 
and (4.10) in the right half-plane. For brevity we write in (4.13) simply A;(cx)) instead of 
lima^oo A;(a2:). Transition to the complex frequencies in eq. (4.13) is accomplished by 
introducing the unphysical photon mass p 

u ^ {z^ + ■ (4.14) 

Extension to the complex 2;-plane of the frequency equations Aj^{az) and A™{az) should 
be done in the following way. In the upper (lower) half-plane the Hankel functions of the hrst 
(second) kind H^\az) {H^\az)) must be used. If we are considering only outgoing waves 
it is natural to use the solutions of the Maxwell equations which are hnite (or at least do not 
grow inhnitely) in the future. Such solutions should be proportional to exp 
when uj lies in the lower half-plane. For the upper half-plane cu, the solutions describing 
outgoing waves and hnite in the future should contain the factor exp 

Location of the roots of eqs. (4.9) and (4.10) enables one to deform the contour C into a 
segment of the imaginary axis (—iA, iA) and a semicircle of radius A in the right half-plane. 
When A tends to inhnity the contribution along the semicircle into Cbau(>s) vanishes because 
the argument of the logarithmic function in the integrand tends in this case to 1. As a result 
we obtain 




^+200 


dz d Aj^{z)A™{z) 


1=1 


2m p^o 


z‘^ + fj?y dz Aj^{ioo)A™{oo) 


(4.15) 


Now we impose the condition that the velocity of light inside and outside the ball is the 
same 
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( 4 . 16 ) 


— £21^2 — c 


The physical motivation for this is the following. The constancy condition for the velocity 
of gluonic held when crossing the interface between two media is used, for example, in a 
dielectric vacuum model (DVM) of quark conhnement [67,125,126]. This model has many 
elements in common with the bag models [127], but among the other differences, in DVM 
there is no explicit condition of the held vanishing outside the bag. It proves to be important 
for calculation of the Casimir energy contribution to the hadronic mass in DVM. The point 
is that in the case of boundaries with nonvanishing curvature there happens a considerable 
(not full, however) mutual cancellation of the divergences from the contributions of internal 
and external (with respect to the boundary) regions. If only the held conhned inside the 
cavity is considered, as in the bag models [128-131], then there is no such a cancellation, and 
one has to remove some divergences by means of renormalization of the phenomenological 
parameter in the model dehning the QCD vacuum energy density. 

/^From a physical point of view the vanishing of the held or its normal derivative precisely 
on the boundary is an unsatisfactory condition, because, due to quantum huctuations, it is 
impossible to measure the held as accurately as desired at a certain point of the space [132]. 

Under condition (4.16) the argument of the logarithm in (4.15) can be simplihed consid¬ 
erably [25] with the result 


Cball('S) 


a\ sinfTTs) 


c/ TT 


^(2/ -h 1) 


1=1 


dy d 

dy 


ln[l (^i (y)], 


(4.17) 


where 


X 2 / \ 2 

-2 / ^1 ~ ^2 \ / hi ~ 7^2 


r = 


, , , My) =-^[si{y) ei{y)]. 

hi + h 2 / dy 


^Si + 62 

Here si{y) and ei{y) are the modihed Riccati-Bessel functions 


(4.18) 


-Si(a^) 



ei{x) 



u = 1 + 1/2. 


(4.19) 


More details concerning the contour integral representation of the spectral ( functions can 
be found in [130,133-135]. 

Further the analytic continuation of eq. (4.17) is accomplished by expressing the sum over 
/ in terms of the Riemann zeta function. This cannot be done in a closed form. Making use 
of the uniform asymptotic expansion (UAE) for the Bessel functions^ in increasing powers 


^In physical literature this expansion is usually referred to as the Debye expansion. However we 
did not find these formulae in Debye papers. Furthermore in mathematical handbook [136] the 
Debye type expansions for the Bessel functions are defined as those for large argument x and large 
order u both real and positive such that v/x is fixed and simultaneously \x — v\ = 0(x^^^). The 
uniform asymptotic expansions were derived in ref. [137] 
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of 1/z/ enables one to construct the analytic continuation looked for in the form of the series, 
the terms of which are expressed through the Riemann zeta function. The problem of the 
convergence of this series does not arise because its terms go down very fast, at least when 
calculating the vacuum energy. 

We demonstrate this keeping only two terms in UAE for the product of the modihed 
Bessel functions Iy{i'z)K,y{i'z) [138] 


Iu{h'z)K^{h>z) ~ 


t 





1 

Vi + 2 :^ 


(4.20) 


After changing the integration variable y = vz in eq. (4.17) we substitute (4.20) into this 
formula and expand the logarithm function up to the order keeping at the same time 
only the terms linear in The last assumption is not crucial. It is introduced for simplicity 
and in order to have possibility of a direct comparison with the results of ref. [25]. Thus we 
have 




-^{zZ{yz)K^{uz)) 


(4.21) 


3 

= -^zt^ + zt^i-12 + 216t2 - 600f^ + 420t®) + C>(z/-®) 

2 z/2 16z/4 ^ \ > 

Integration over 2 ; can be done by making use of the formula [112] 


A(^')"(-f[ 

2 : “ ^t^dz = -—-, Rea < 0, Re (a + /?) > 0. 

T f ^ 

2 


(4.22) 


Also the properties of the T function 


r(;.)r(i-;.) = 


TT 


r(l + ;^) = ;2r(;2) 


SlUTT^: 

prove to be useful. After simple calculations we arrive at the result 


(4.23) 


Cb.ii(s) ^ (((-)'%(!+ s )(2 + s) I 5 ^ + p(s) Y, + ■■■). (4.24) 


. 1=1 


1=1 


where 


P(^) = 


1 - -(3 + s) + -(3 + s)(4 + s) - —(3 + s)(4 + s)(5 + s) 


(4.25) 


The hrst term in eq. (4.24) is dehned in the region 0 < Res < 1/2 and the second one 
exists, when —1 < Res < 1/2, the left boundaries being given by the sums over the values 
of the angular momentum I and the right limits are due to the integral formula (4.22). Thus, 
this expression can be used for its analytical continuation outside these regions. Obviously, it 
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can be done by accepting the integration formula (4.22) outside the domain indicated there 
and by expressing the sums over angular momentum I trough the Riemann zeta function 
according to the formula [ 112 ] 


z/-* = (2^ - 1)Cr(s) - 2% u = 1 + 1/2. (4.26) 

1=1 

As a result one gets 

Cbaii(s) ~ ^ s(l + s)(2 + s){(2i+2* _ + 2s) - 2 ^+ 2 . 

+p{s) [(2^+2^ - 1)Cr(1 + 2s) - 2^+2^] + ...}. (4.27) 

The singularities in eq. (4.24) are transformed in (4.27) into the poles of the Riemann zeta 
functions at the points s = —k, k = 0,1,2 ,... 

Cr(1 + 2s)~^ + 7 + ..., s^O, 

Cr(3 + 2s) ~ ^^+7 + ..., s^-2, (4.28) 

Cr(5 + 2s) ~ 23^+4 + 7' + • • • > 8^-4, 


where 7 is the Euler constant. The first three poles are annihilated by the multipliers in front 
of the curly brackets in eq. (4.27). The nearest surviving singularity (simple pole) appears 
only at the point s = —3. Thus formula (4.27) affords the required analytic continuation of 
the function CbaiK-s) into the region Re s < 0. In view of eq. (3.13) we are interested in the 
point s = —1/2 where CbaiK-s) given by (4.27) is regular 


A-ball — 



64a 


1 + 


9 


T-4 

128 V 2 


+ 


(4.29) 


It is exactly the hrst two terms in eq. (3.10) of ref. [25]. The procedure of analytic continua¬ 
tion presented above can be extended in a straightforward way to the arbitrary order of the 
uniform asymptotic expansion (4.20). Certainly in this case analytical calculations should 
be done by making use of Mathematica or Maple. 

In sect. V the exact in the .^^-approximation expression for the Casimir energy (4.29) will 
be obtained by making use of another approach (see eq. (5.27)). 

The problem under consideration with .^ = 1 is of a special interest because in this 
case it gives the Casimir energy of a perfectly conducting spherical shell. As it was noted 
above, this conhguration has been considered by many authors. We present here the basic 
formulae which afford the analytical continuation of the corresponding spectral zeta function. 
We again content ourselves with two terms in the UAE (4.20). It is impossible to put simply 
(^ = 1 in the next formula (4.21). One has to do the expansion here anew keeping all the 
terms ~ 1/z/"^. This gives 
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1 - 



{zIi,{uz)K^{uz)) 



= ^zf + -^zt^ (-1 + 18f - 50t^ + 35t®) + C>(z/-®) . (4.30) 

2z/^ 4z/^ J 

After integration and elementary simplifications we arrive at the following result for the 
spectral function in hand 

2 s oo oo 

Csheii(s) ^ + s) (2 + s) + q{s) Y + • • • > (4-31) 

li=i 1=1 


where 

q(s) = ^ (60 + 217s + 2525^ + 71s^). (4.32) 

^ ^ 480 

The analytic continuation of eq. (4.31) is accomplished in the same way as it has been done 
with eq. (4.24), i.e., , by making use of eqs. (4.22) and (4.26) 

CshelK-S) — > 5(1 + s)(2 + s){(2^’''^^ — 1 )Cr( 1 + 2s) — 2^^^® + 

+g(s) [(2^+2^ - 1 )Cr( 3 + 2s) - 2^+2®] + ...}. (4.33) 


The nearest singularity in this formula is simple pole at s = —3. As above it is originated in 
the term ~ 1/z/’^ in the UAE (4.20). At the point s = —1/2 the spectral zeta function Csheii('S) 
is regular and gives the following value for the Casimir energy of a perfectly conducting 
spherical shell 


^sheii 2 ) 64a ^ 256(2 A+--- 


= -0.046361 
a 


(4.34) 


Without considering the analytic continuation and do not carrying out the analysis of the 
singularities in the complex s plane this result has been obtained in [118]. For nanometer 
size, that is for a = 10“'^cm, the energy (4.34) is (in h = c = 1 unit, leV ~ 0.5 ■ lO^crn”^) 
-E'sheii — lOeV which is of a considerable magnitude. 

Unlike the Casimir energy for parallel conducting plates (4.5) both the Casimir energies 
for dielectric-diamagnetic ball (4.29) and for perfectly conducting spherical shell (4.34) prove 
to be positive. It implies that the Casimir forces in these cases are repulsive. In view of 
this the Casimir idea [139] to construct an extended model for electron failed because the 
vacuum forces cannot stabilize it. 

In principle, it is possible to get a negative Casimir energy for a sphere, however for this 
purpose we have to consider the scalar massless field obeying on inner and outer surfaces of 
the sphere the Neumann boundary conditions [118] 

= --0.223777.... (4.35) 

a 
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The same field obeying the Dirichlet conditions on a sphere has positive vacuum energy [118] 

=-0.002819.... (4.36) 

a 

When obtaining the Casimir energies (4.35) and (4.36) the summation over the values of the 
angular momentum I starts from / = 0 , unlike the electromagnetic field, where I assumes 
the values 1 , 2 ,.... 

Undoubtedly, the calculation of the Casimir energy for a material ball with ei/ii 7 ^ £ 2/^2 
by a rigorous zeta function method is also of a special interest. However, in this case the very 
definition of the spectral zeta function should probably be changed in order to incorporate 
the contact terms which seem to be essential in this problem [38,40,41,140] (see sect. V). 


C. Cylinder 

Calculation of the Casimir energy of an circular infinite cylinder [20,19] proves to be a 
more involved problem in comparison with that for sphere. Here the spectral zeta function 
CcyK-s) for this configuration will be constructed, its analytical continuation into the left 
half-plane of the complex variable s will be carried out, and relevant singularities will be 
analyzed. 

Thus we are considering an infinite cylinder of radius a which is placed in an uniform 
unbounded medium. The permittivity and the permeability of the material making up 
the cylinder are ei and pi, respectively, and those for surrounding medium are 62 and fi 2 - 

We assume again that the condition (4.16) is fulfilled. In this case the electromagnetic 

oscillations can be divided into the transverse-electric (TE) modes and transverse-magnetic 
(TM) modes. In terms of the cylindrical coordinates (r, 6 , z) the eigenfunctions of the given 
boundary value problem contain the multiplier 

exp {±iut + ikzZ + inO) (4.37) 

and their dependence on r is described by the Bessel functions Jn for r < a and by the Hankel 
functions or for r > a. The frequencies of TE- and TM-modes are determined, 
respectively, by the equations [ 122 ] 

A™(A, a) = Xa[niJ^{Xa)Hn{Xa) — fiiJn{Xa)H'^{Xa)] = 0 , (4.38) 

A™(A, a) = Xa[eiJ^{Xa)Hn{Xa) - eMXa)H^{Xa)] = 0 , (4.39) 

where n = 0,±1,±2,.... Here A^ is the eigenvalue of the corresponding transverse 
(membrane-like) boundary value problem 

X^ = ^-kl (4.40) 

Decoupling of TE- and TM-modes in this problem is due to the condition (4.16). 

In a complete analogy with the ball the Casimir energy per unit length of the cylinder is 
defined by the relevant spectral zeta function according to eq. (3.13). Let Xnr be the roots 
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of the frequency equations (4.38) and (4.39), then the function Ccyi(’S) is introduced in the 
following way 


r*+oo 


CcyK-S) — C 


-2s 


dkz 


- (y,r(oo) + *1 


(4.41) 


In terms of the contour integral (see eq. (3.16)) it can be represented in the form 

r A™(Aa)A™(Aa) 


Ccyl(s) = C ^ ^ 

J—oo C” 


ATE(oo)A™(c 5 o) ■ 


(4.42) 


Again we can take the contour C to consist of the imaginary axis (+fcxD, —icx)) closed by a 
semicircle of an inhnitely large radius in the right half-plane. Continnation of the expressions 
A™(Aa) and A™(Aa) into the complex plane A should be done in the same way as in the 
preceding section, i.e., by using Hn\X) for Im A < 0 and Hn\X) for Im A > 0. On the 
semicircle the argument of the logarithm in eq. (4.42) tends to 1. As a resnlt this part of 
the contour C does not give any contribntion to the zeta fnnction Ccyi('S). When integrating 
along the imaginary axis we choose the branch line of the fnnction /(A) = (A^ + A;^)”^ fo rnn 
between —ik^ and +ikz, where kz = +a/^ > 0 and use further that branch of this function 
which assnmes real values when \y\ < kz, where y = Im A. More precisely we have 

( e-“®(|/2 - kl)-% y > kz, 

f{iy) = I {kl-y‘^)-% \y\ < kz, (4.43) 

[ - kl)-\ y < -kz. 

Employment of the Hankel functions Hn\\) and Hn\X) by extending the expressions 
A™(A) and A™(A) into the complex plane A, as it was noted above, gives rise to the 
argument of the logarithm fnnction depending only on on the imaginary axis. It means 
that the derivative of the logarithm is odd fnnction of y. As a result the segment of the 
imaginary axis {—ikz, +ikz) gives zero, and eq. (4.42) acqnires the form 


Ccyi(s) 


sm ns 


q2s^2 


+ 00 
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n=—oo 




k‘1) ^dyhi 


AT^^(a|/)„A™(a|/) 

A™(ioo)A™(ic5o) ■ 


(4.44) 


Changing the order of integration of kz and y and taking into acconnt the value of the 
integral [112] 


i_2s r(i-s) 


/ dkz{y^-kir = ^y 

^ r --S 


Re s < 1, 


(4.45) 


we obtain after the snbstitution ay ^ y 


CcyK-S) 


2v/7rar(s)r ( 2 “ ^ 


2 ”l“00 ^ oo ^ 

0 I In[1(?/)]> (4-46) 
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where 


^^n{y) = y[In{y)Kn{y)\, (4.47) 

and the parameter was defined in eq. (4.18). We shall again content ourselves with 
the first two terms in the uniform asymptotic expansion (4.20) and take into account only 
the terms linear in In this approximation, upon changing the integration variable y = 
nz, n = ±1, ±2 ,..., we have 


1 

f 2 

■ d 

2] 

In < 


z—{In{nz)Kn{nz)) 

dz 

J 


(4.48) 


= -e‘ 


zH^ 

4n^ 


1 + ^(3 - 30r + 35r) + 0{n-*) 


Now we substitute (4.48) into all the terms in (4.46) with n ^ 0. The term with n = 0 in 
this sum will be treated by subtracting and adding to the logarithmic function the quantity 


y^ 


4 (1 + 2/2)3 ■ 

As a result the zeta function CcyK-s) can be presented now as the sum of three terms 


CcyK"®) •^l('®) T •^2('S) + ) 


(4.49) 


(4.50) 


where 


ZAs) = 


2s 


2y/¥ar(s)r ( 2 “ ^ 


7 “ d r 72 


(4.51) 


+ CXD 


^ 2 ( 5 ) = (- 


2s 


2 ^n“ 2 *-i + 1 


n=l 


8 v^ar(s)r 

! E 
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/ dzz^-^ — {zH^), 
I Q dz 


(4.52) 


-3-2s 


^3(s) = 


n=l 


00 . 

l-2s “ 


?>2Z^a r(s)r 


dzz^-‘^^^[zH^{‘i - 3072 + 357 ^)] . 
^0 dz 


(4.53) 


In these equations Zi(s) has accumulated the term with n = 0 from eq. (4.46) subtracted 
by (4.49); ^ 2 ( 5 ) involves the contribution of the term of order l/n2 in expansion (4.48) and 
the added expression (4.49); ^ 3 ( 5 ) is generated by the terms of order l/n^ in the expansion 
(4.48). 

Taking into account that 
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/^o(l/)ly-o ^ 1 and ^il{y)\y^^ ^ — + 


(4.54) 


4|/2 IQy* ’ 

the integration by parts in eq. (4.51) can be done for —3/2 < Re s < 1/2 with the result 
(2s - l)a2*-i 


ZAs) = 


2Virc2*r(s)r 




/o y 


(4.55) 


With allowance for (4.54) one infers easily that the function Zi{s) is an analytic function of 
the complex variable s in the region —3/2 < Re s < 1/2. In the linear order of it reduces 
to 


Zl-(s) = f 


2s - 1 


2v/iar(s)r ( 2 “ 


a\ 


2s 


dy 


c/ Jo y 


.2s 
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4(1 + 1/2)3 


- yliy) 


(4.56) 


This function is also analytic in the region —3/2 < Re s < 1/2. Integration in eq. (4.52) 
can be accomplished exactly by making use of the formula [112] 


dz z 
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l-2s /„4.6 


dz 


z‘^n = 


2s- 1„ fl 


+ ,-J<Res<^. (4.57) 


This gives for Z 2 {s) in (4.52) 


z^{s) = e(-Y^ 
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64^/na 
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-2s-l 


+ 1 


n=l 


ri2+* 

r(i>) 


(4.58) 


In view of the sum over n in (4.58) the function Z 2 {s) is dehned only for s > 0. 

For simplicity we apply in eq. (4.53) the integration by parts which is correct for —3/2 < 
Re s < 1 and leads to the result 


rf^o\2*(l-2s)(3-2s)(28s2 

hi 6144v4Fa 


8s - 27) 


+ s 


r(s) 


n=l 


- 25-3 


(4.59) 


Again the sum over n in (4.59) gives the restriction Re s > —1 for the dehnition of the 
function Z 3 (s). 

Thus the spectral zeta function Ccy^s) in the linear approximation with respect to and 
with allowance for the first two terms in the UAE (4.48) is given by 

Ccyi(s) = Zf(s) + Z2{s) + Z3(s), (4.60) 

where the Z’s are presented in eqs. (4.56), (4.58) and (4.59), respectively. Summing up 
all the restrictions on the complex variable s which have been imposed when deriving these 
equations we infer that CcyK-s) is dehned in the strip 0 < Re s < 1/2. In order to continue 
this function into the surroundings of the point s = —1/2, it is sufficient to express the 
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sum in eq. (4.59) in terms of the Riemann zeta function and treat the right-hand side of eq. 
(4.57) as an analytic continuation of its left-hand side over all the complex plane s. This 
gives 


Z,(s) 


^2 ("“V" (1 -2s)(3- 2s) 
\cJ 6A^/^T a 


r 

[2Cr(2"S + 1) + 1] — 



(4.61) 


It is left now to take the limit s —1/2 in eqs. (4.56), (4.58) and (4.61). A special care 
should be paid when calculating this limit in eq. (4.61) in view of the poles of the function 
r(l/2 -|- s) at this point. Using the formulae 

Cr(0) =Cr(0) =-^ln27r, T{x) = - -'-f + 0{x) a; ^ 0 , (4.62) 

2 2 X 

one derives 


lim [2 Cr( 1 -h 2s) -h l]r 

s^-l/2 


1 + 2s 


= lim [2Cr(0)+2C(,(0)(1 + 2s) + 0((1 + 2s)2) + 1] 

S —^ 1 ^ 

Making allowance for this, we obtain from (4.61) 


1 + 2s 


+ 7 + 0(1 +2s) 


Z2(-l) 
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(4.63) 


- 21 n( 27 r). 


(4.64) 


The appearance of the hnite term proportional to ln(27r) is remarkable for the problem 
under consideration. It is derived here in a consistent way by making use of an analytic 
continuation of the relevant spectral zeta function. In ref. [19] it was calculated in a more 
transparent though not rigorous way. 

Gathering together eqs. (4.56), (4.59) with s = —1/2 and eq. (4.64) we obtain 


Cc,i(-l/2) = 


+1 

27ra^ 


ydy 
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[A{l + y^y 


-yliy) 


1^1 1,2, 

+ — ^ —5- + — ln( 27 r) 

48 ^ 4 

n=l 


= (-0.490878 + 0.034269 + 0.459469):/^ = 0.002860 ^ 


27ra^ 


27ra^ 


(4.65) 


This result is not the hnal answer in the problem in hand. Point is that in view of severe 
cancellations in eq. (4.65) the contribution of the next term in the UAE (4.48) proves to be 
essential. Its account gives [19] 


Ccyi(-l/2) = 0. (4.66) 

Thus the Casimir energy of a compact cylinder possessing the same speed of light inside and 
outside proves to be zero in the ^^-approximation. The consideration presented in this section 
can be extended to the next term of order ~ 1/n® in the UAE (4.48) in a straightforward way. 
Therefore we shall not present here these rather cumbersome expressions. More precisely, in 
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ref. [19] not only the next term in UAE (4.48) has been taken into account but also the first 
5 terms in the sum in (4.53) were taken exactly instead of using their asymptotics . 

This turns out to be essential for reaching the required precision. However this point is not 
critical for our consideration because Z^ls) does not need analytic continuation. The exact 
zero result for this energy was obtained in paper [35] by making use of the addition formulae 
for the Bessel functions instead of the UAE. 

However the Casimir energy of a compact cylinder possessing the same speed of light 
inside and outside does not vanish in higher approximations with respect to the parameter 
[26,35]. For example, up to the term it is dehned by the formula 

E{^^) = -0.0955275^^ = -0.007602 (4.67) 

4t7r Qj Qj 

In contrast to the Casimir energy of a compact ball with the same properties (see eq.(4.29)) 


A-ball — 


= 0.046875 — 
64a a 


(4.68) 


the Casimir energy of a cylinder under consideration turned out to be negative. Conse¬ 
quentially, the Casimir forces strive to contract the cylinder. The numerical coefficient in 
eq. (4.67) proved really to be small, for example, in comparison with the analogous coeffi¬ 
cient in eq. (4.68). Probably it is a manifestation of the vanishing of the Casimir energy of 
a pure dielectric cylinder noted in the Introduction. Numerically the Casimir energy 
was calculated for arbitrary values of the parameter in paper [26]. 

Now we address to the consideration of a special case when = 1. It corresponds to a 
perfectly conducting cylindrical shell [19]. Instead of the expansion (4.48) we have 
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1 
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(4.69) 


4n^ 


1 + ^(3- 30C + 35t^ + - ) + 0(71-^) 


Proceeding in the same way as above we obtain for the spectral zeta function concerned 


CS;,"'(s) = Z,(s) + Z2(s) + Z:t(3) , (4.70) 

where Zi{s) is given by eq. (4.55) with = 1; Z 2 (s) is the same as in eq. (4.61), and Z^ls) 
now is 


^ 3 ( 5 ) 
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+ 00 
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(4.71) 


At the point s = —1/2 it acquires the value 
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87ra^ 
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= -0.0268- 




(4.72) 



This exactly reproduce the contribution of the first two terms in calculations of the Casimir 
energy for perfectly conducting circular cylindrical shell in ref. [19]. With higher accuracy 
this energy is given by [19,20] 

R^ell ^ _Q_Ql35gi__ 

The Casimir energy of a massless field subjected to Dirichlet and Neumann conditions 
on internal and external surfaces of an inhnite circular cylindrical shell are also hnite but of 
opposite signs [141] 

= 0 . 0006064 , 

CL 

El, = -0.014174- (4.75) 

CL 

The sum of these two energies gives the vacuum energy of electromagnetic held for this 
conhguration (4.73). 

However the zeta regularization does not give a hnite value for the Casimir energy in 
the case of the two-dimensional version of the conhguration under consideration, i.e., for a 
circle of radius a placed on a plane. In paper [141] these energies were calculated by making 
use of the relation between the zeta functions for a cylinder and a circle of the same radius 

Ccir(s) = + 1/2). (4.76) 

Having calculated the zeta function for an inhnite cylinder in the region —1/2 < Res < 0 
one obtaines immedeatly the Casimir energy for a cylinder via /’cyi(—1/2) and the Casimir 
energy for a circle in terms of Ccyi(O) 

El = - (0.0023595 - —- 
a V 256 s 


(4.77) 
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The Casimir energy for the sum of these two helds is also inhnite 


(4.78) 
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(4.79) 


On a plane electromagnetic held is reduced to a scalar massless held obeying the Neumann 
boundary conditions [142]. It was shown, that the zeta function regularization does not 
lead to a hnite answer for the vacuum energy for spheres in spaces of arbitrary even dimen¬ 
sions [142-145]. As usual the coefficients in front of the pole-like contributions calculated 
by diherent methods coincide, but the hnite parts of the answers diher [142,143]. Thus for 
obtaining here physical results additional renormalization is needed. 
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In ref. [21] the vacuum energy of a perfectly conducting cylindrical surface has been 
calculated to much higher accuracy by making use of another version of the zeta function 
technique. By integrating over dk^ directly in eq. (4.41) the authors reduced this problem to 
investigation of the zeta function for circle, which has been considered earlier by introducing 
the partial wave zeta functions for interior and exterior region separately. In this respect our 
approach dealing only with one spectral zeta function for given boundary conditions proves 
to be more simple and straightforward. 

Many results on calculation of the Casimir energies for spherically symmetric cavities by 
making use of the zeta function technique can be found in ref. [145]. 


V. CASIMIR ENERGY OF A DIELECTRIC BALL 

Calculation of the vacuum energy of electromagnetic held connected with a dielectric ball 
proves to be a quite involved problem. The hrst attempt to calculate this energy has been 
undertaken by K. A. Milton in 1980 [37]. And only just recently the vacuum electromagnetic 
energy of a dilute dielectric ball was found at zero [38-45] and hnite [46,47] temperature. It 
is worth noting that the hrst rough and not rigorous estimation of the Casimir energy of a 
dilute dielectric ball has been done in ref. [25]. The difficulties encountered here are due to 
the following mathematical peculiarities of the problem in question. 

As it was noted in subsect. IV B the solutions to the Maxwell equations with spherical 
symmetry are expressed in terms of two scalar Debye potentials obeying Helmholtz equa¬ 
tions (4.8). In the Casimir calculations these equations should be treated as one eigenvalue 
problem, hence we have to rewrite them in the form 

c2 

i = 1 for r < a, i = 2 for r > a. (5.1) 

Here Ci and C 2 are the velocities of light respectively inside and outside the ball. Thus the 
coefficient in front of the differential operator —A is discontinuous, it has a jump at the 
surface of the ball. It turns out that the heat kernel coefficient B 2 for the spectral problem 
(5.1) does not vanish but is a function of ci — C 2 [146]. If this difference is small (a dilute 
dielectric ball) it was shown that 


B 2 ~ (ci - C2)^ (5.2) 

Thus the coefficient B 2 vanishes up to terms (ci — 02 )^. It is this approximation that have 
been used in all the calculations of the Casimir energy of a dielectric ball. 

In sect. HI we have tried to show that the zeta function technique is the most consistent, 
from the mathematical standpoint, method for calculation of vacuum energy. However one 
did not manage to use this technique in a direct way for hnding the Casimir energy of a 
dielectric ball. Here we describe briefly the methods practically used in these calculations. 

In ref. [38] the Casimir energy of a dilute dielectric ball was derived by summing up the 
van den Wails forces between the individual molecules inside the compact ball. Barton [39] 
found this energy in the framework of a special perturbation theory, in which the dielectric 
ball was treated as a perturbation when treating the electromagnetic held in unbounded 
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empty space. Milton, Brevik and Marachevsky [40,41,43] started from the Green’s function 
of the quantized Maxwell field with an explicit account of the so called contact terms, on the 
stage of the numerical calculations the uniform asymptotic expansion for the Bessel functions 
and the zeta regularization technique being applied. We have considered this problem [42] by 
direct mode summation method with use of the addition theorem for the Bessel functions. 
In calculations without employment of the uniform asymptotic expansions for the Bessel 
functions [38,39,42] the exact (in the An^ approximation) result for the Casimir energy at 
hand was obtained. Brevik et al [44,45] considered this problem in a statistical approach. 

First we present the calculation of this Casimir energy at zero temperature and further 
we extended this approach to a hnite temperature. 


A. Zero temperature 

As in subsect. IV B, a solid ball of a radius a placed in an unbounded uniform medium 
is considered but the condition (4.16) is not imposed now. Unfortunately, the zeta regular¬ 
ization cannot be used here in a direct way. The reason for this will be seen below. 

We shall proceed from the standard dehnition of the vacuum energy as the sum over the 
eigenfrequencies of electromagnetic oscillations 

E = ^'^iuJs-Us). (5.3) 

Here Ug are the frequencies dehned by eqs. (4.9) and (4.10), and the frequencies Ug correspond 
to a certain limiting boundary conditions that will be specihed below. 

In the case of the plane geometry of boundaries or when considering the Casimir effect 
for distinct bodies it is sufficient to subtract in eq. (5.3) the contribution of the Minkowski 
space. In the problem at hand it implies taking the limit a oo, i.e., that the medium 1 
tends to £11 the entire space. But it turns out that this subtraction is not sufficient because 
the linear in ei — 62 contribution into the vacuum energy retains. Further, we assume that 
the difference ei — 82 is small and content ourselves only with the (si — e 2 )^-terins. 

The necessity to subtract the contributions to the vacuum energy linear in ei — 82 is 
justified by the following consideration. The Casimir energy of a dilute dielectric ball can be 
thought of as the net result of the van der Waals interactions between the molecules making 
up the ball [38]. These interactions are proportional to the dipole momenta of the molecules, 
i.e., to the quantity (ei — 1)^. Thus, when a dilute dielectric ball is placed in the vacuum, 
then its Casimir energy should be proportional to (si — 1)^. It is natural to assume that 
when such a dielectric ball is surrounded by an infinite dielectric medium with permittivity 
82 , then its Casimir energy should be proportional to (ei — 82 )“^- The physical content of the 
contribution into the vacuum energy linear in 81 —82 has been investigated in the framework 
of the microscopic model of the dielectric media (see paper [147] and references therein). It 
has been shown that this term represents the self-energy of the electromagnetic held attached 
to the polarizable particles or, in more detail, it is just the sum of the individual atomic Lamb 
shifts. Certainly this term in the vacuum energy should be disregarded when calculating 
the Casimir energy which is originated in the electromagnetic interaction between different 
polarizable particles or atoms [39,40,44,45,148]. 
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Further, we put for sake of symmetry 


^ = n 2 = 1 


(5.4) 


^ = ni 



An 

~Y 


Here ni and n 2 are the refractive indices of the ball and of its surroundings, respectively, 
and it is assumed that An << 1. From here it follows, in particular, that 


£i- £2 = + n2)(ni - n2) = 2An. 


(5.5) 


Thus, using the dehnition (5.3) we shall keep in mind that really two subtractions should 
be done: hrst the contribution, obtained in the limit a 00 , has to be subtracted and then 
all the terms linear in An should also be removed. 

We present the vacuum energy dehned by eq. (5.3) in terms of the contour integral in 
the complex frequency plane. Upon the contour deformation one gets 
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^( 2 / + 1 ) 

1=1 



dyy — \n 
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A™ {iay)A™ {iay) 
A™(icx))A™(icx)) ’ 


(5.6) 


where A™(my) and A™(my) are dehned in eqs. (4.9) and (4.10). 

In eq. (5.6) we have introduced cutoff K in integration over the frequencies. This regu¬ 
larization is natural in the Casimir problem because physically it is clear that the photons of 
a very short wave length do not contribute to the vacuum energy since they do not “feel” the 
boundary between the media with different permittivities ei and 62 - In the hnal expression 
the regularization parameter K should be put to tend to inhnity, the divergencies, that may 
appear here, being cancelled by appropriate counter terms. 

Taking into account the asymptotics of the Riccati-Bessel functions we obtain 

A™(, 00 )A™(^ 00 ) = _ e 2 (ni-n.), _ ^5 7 ) 

Upon snbstitnting eqs. (4.9), (4.10) and (5.7) into eq. (5.6) and changing the integration 
variable ay y, we cast eq. (5.6) into the form (see ref. [25]) 
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(5.8) 


where si = siijiiy), ei = ez(n 2 y), yo = aK. 

It shonld be noted here that in eq. (5.8) only the first snbtraction is accomplished, which 
removes the contribntion to the vacunm energy obtained when a ^ oo. As noted above, for 
obtaining the hnal resnlt all the terms linear in An should be discarded also. 

Further it is convenient to rewrite eq. (5.8) in the form 


E — El + E 2 


(5.9) 


with 
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( 5 . 10 ) 
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ydy, 
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(5.11) 


(5.12) 

(5.13) 


The term Ei accounts for only the expression exp (—2An y) in the argument of the loga¬ 
rithm function in eq. (5.8) and it appears as a result of subtracting the Minkowski space 
contribution to the Casimir energy (the sum with ujg in eq. (5.3) and the denominator in eq. 
(5.6)). 

It is worth noting that the term Ei is exactly the Casimir energy considered by Schwinger 
in his attempt to explain the sonoluminescence [149]. Really, introducing the cutoff K for 
frequency integration and the cutoff y = uj/a for the angular momentum summation we 
arrive at the result 
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(5.14) 


We have substituted here the summation over / by integration. Up to the multiplier (—2/3) 
it is exactly the Schwinger value for the Casimir energy of a ball (ei = 1) in water (y^ — 
4/3) [140]. The term linear in An and of the same structure was also derived in papers 
[39,44,148]. As it was explained above the energy Ei should be discarded. 

In our calculation, we content ourselves with the An^-approximation. Hence, in eq. 
(5.11) one can put P^{niy,n 2 y) — P^{y,y) and keep in expansion of the logarithm function 
only the terms proportional to An^. In this approximation, the contributions of and P/ 
into the vacuum energy are additive 


^ren _ ^ 

In the Appendix A it is shown that for obtaining the An^-contribution into the Casimir 
energy of the function W/ in the argument of the logarithm in eq. (5.11), it is sufficient 
to calculate the An^-contribution of the function alone but changing the sign of this 
contribution to the opposite one (see eq. (A20)). Hence, 

1 d 

+ 1) J dyy—Wf{niy,n2y), (5.16) 


and only the An^-term being preserved in this expression. 

For Ep we have 

A972 ryo y 

Ep = -— ^(2/ + 1) / dyy—Pf(niy,n 2 y) 

STia ^ Jq dy 


(5.17) 
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Usually, when calculating the vacuum energy in the problem with spherical symmetry, the 
uniform asymptotic expansion of the Bessel functions is used (see sect. IV). As a result, an 
approximate value of the Casimir energy can be derived, the accuracy of which depends on 
the number of terms preserved in the asymptotic expansion. 

We shall persist in another way employing the technique of the paper [33]. By making 
use of the addition theorem for the Bessel functions [138], we hrst do the summation over the 
angular momentum I in eq. (5.11) and only after that we will integrate over the imaginary 
frequency y. As a result, we obtain an exact (in the An^-approximation) value of the Casimir 
energy in the problem involved. 

The addition theorem for the Bessel functions is given by the formula [138] 

OO ^ 

y^(2/ + l)si{Xr)ei{Xp)Pi{cose) = (5-18) 

1=0 


where 


R = A/r^ + _ 2rp cos 9 . 


(5.19) 


Differentiating the both sides of eq. (5.18) with respect to Xr and squaring the result we 
deduce 
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Here the orthogonality relation for the Legendre polynomials 

J ^ Pl{x)Pm{x)dx = 

has been taken into account. Now we put 

, , An An 

X = y, r = ni = l + —, p = n 2 = l- —. 

Applying eq. (5.20) and analogous ones, we derive 
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(5.22) 

(5.23) 


Here and Pp stand for the results of the partial differentiation of the function P in eq. 
(5.18) with respect to the corresponding variables and with the subsequent substitution of 
(5.21). The last terms in eqs. (5.22) and (5.23) are WQ{niy,n 2 y) and P^iy^y), respectively. 
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As it was stipulated before, in eq. (5.22) we have to keep only the terms proportional to 
An^ and in eq. (5.23) we have put An = 0. 

The calculation of the contribution Ep to the Casimir energy is straightforward. Upon 
differentiation of the right-hand side of eq. (5.23) with respect to y, the integration over dR 
can be done here. Substitution of this result into eq. (5.17) gives 


Ep = 





+ 2 |/ + 


1 
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(5.24) 


The term (—1/2) in the square brackets in eq. (5.24) gives rise to the divergence^ when 
I/O ^ oo 
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div 
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Ibvra 
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(5.25) 


Therefore we have to subtract it with the result 

^ /\ti^ 

Er = Ep-Ef^ = — -. (5.26) 

12o TTO 

It is worth noting here that eq. (5.26) after substitution An^/4 = gives the exact, in 
the .^^-approximation, value for the Casimir energy of a compact ball with the same velocity 
of light inside and outside the ball [33,42,46,15] 


-E-ball — 
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(5.27) 


(compare this formula with eq. (4.29)). In subsect. IVB this energy has been calculated 
by making use of the zeta function technique (see eq. (4.29) ). As this zeta function is not 
known exactly, we have derived there only an approximation for the exact result (5.51). 

As far as the expression (5.22), it is convenient to substitute it into eq. (5.16), to do the 
integration over y and only after that to address the integration over dR 
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(5.28) 


We have put here yo = oo without getting the divergencies. As it is explained in the 
Appendix A, in eq. (5.28) we have to pick up only the term proportional to An^. Remarkably 


^This divergence has the same origin as those arising in summation over I when the uniform 
asymptotic expansions of the Bessel functions are used [40,41]. The technique employed here is 
close to the multiple scattering expansion [15], where these divergencies are also subtracted. 
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this term is finite. It is an essential advantage of our approach. The rest of the terms in this 
equation are irrelevant to our consideration. Thus the counter term for vanishes due to 
the regularizations employed (see the Appendix A). In view of this we have 
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1 1 An^ 1 An^ 


w 


27ra 8 3 48 vra 

Finally we arrive at the following result for the Casimir energy of a dilute dielectric ball 


(5.29) 


^ren 


Tpren , jpre 

-^W ' 


An^ / I 5 \ 

V48 ^ my 


23 An^ 
384 Tra 


(5.30) 


Taking into account the relation (5.5) between and i = 1, 2, we can write 


23 (si — £ 2 )^ 
1536 Tra 


(5.31) 


At the first time, this value for the Casimir energy of a dilute dielectric ball has been de¬ 
rived in ref. [38] by summing up the van der Waals interactions between individual molecules 
making up the ball (^2 = !)• The result (5.31) was obtained also by treating a dilute di¬ 
electric ball as a perturbation in the complete Hamiltonian of the electromagnetic field for 
relevant configuration [39]. In papers [40,41], the value close to the exact one, eq. (5.31), 
has been obtained by employing the uniform asymptotic expansion of the Bessel functions. 

In ref. [25] the estimation of the Casimir energy of a dilute dielectric ball has been done 
taking into account, as it is clear now, only the second term in eq. (5.30). And nevertheless 
it was not so bad having the accuracy about 35%. 


B. Finite temperature 

An essential advantage of the calculation of the Casimir energy of a dilute dielectric ball, 
carried out above by the mode summation method, is the possibility for its straightforward 
generalization to the finite temperature. The employment of the addition theorem for the 
Bessel functions again enables one to carry out the summation over the angular momentum 
in a closed form. As a result, the exact (in the An^ approximation) value for internal 
and free energies, as well as for entropy of a dilute dielectric ball are derived for finite 
temperature also. The divergencies, inevitable in such studies, are removed by making use 
of the renormalization procedure developed for calculation of the relevant Casimir energy 
at zero temperature (see subsect. VA). The thermodynamic characteristics are presented 
as the sum of the respective quantity for a compact ball with uniform velocity of light and 
an additional term which is specific only for a pure dielectric ball. The behavior of the 
thermodynamic characteristics in the low and high temperature limits is investigated. 

Practically the extension to finite temperature T is accomplished by substituting the y- 
integration in eqs. (5.11) by summation over the Matsubara frequencies ujn = 2'KnT. Doing 
in this way we obtain the internal energy of a dielectric ball 
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U(T) = -TJ2(21 + 1)J2' 

1=1 n=0 
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where 


Wi{niWn,n2Wn) = Si{niWn)e[{n2Wn) - s[{niWn)ei{n2Wn) , 
Pl{niWn,n2Wn) = Si{niWn)e\{n2Wn) + s\{niWn)ei{n2Wn) , 

and we have introduced the dimensionless Matsubara frequencies 

Wn = aUn = 2nnaT, n = 0,1, 2,... . 


(5.32) 


(5.33) 

(5.34) 


(5.35) 


The prime on the summation sign in eq. (5.32) means that the n = 0 term is counted with 
half weight. In what follows it turns out to be important, for example, when removing the 
divergencies, that we proceed from the expression for the internal energy instead of from 
the analogous formula for free energy. 

When considering the low temperature behavior of the thermodynamic functions of a 
dielectric ball the term proportional to in our paper [46] was lost. It was due to the 
following. We have introduced the summation over the Matsubara frequencies in eq. (3.20) 
under the sign of the /^-integral. Here we show how to do this summation in a correct way. 
In the A^-approximation the last term in eq. (3.20) from the article [46] 

_ °° r2 „-2w„R 

Uw{T) = 2TAn^2^'^n / —^ Wn = 2nnaT (5.36) 

n=o ^ 

can be represented in the following form 

OO 

UwiT) = -2TAn^ Y^wl Si(4w;„), (5.37) 

n=0 

where Ei{x) is the exponential-integral function [138]. Now we accomplish the summation 
over the Matsubara frequencies by making use of the Abel-Plana formula 

= r f(x) dx + ,r /(“> - /(-»") ,5.38) 

^ J« Jo e"” - 1 

The hrst term in the right-hand side of this equation gives the contribution independent of 
the temperature, and the net temperature dependence is produced by the second term in 
this formula. Being interested in the low temperature behavior of the internal energy we 
substitute into the second term in eq. (5.38) the following expansion of the function Ei{z) 

OO / "I N ^ 

Ei{z) = -7 - In^ - ^ fc . fcf ’ (^-39) 
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where 7 is the Euler constant [138]. The contribution proportional to is produced by 
the logarithmic term in the expansion (5.39). The higher powers of T are generated by the 
respective terms in the sum over k in this formula {t = 27raT) 
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(5.40) 


All these terms, safe for 2(^(3)An^a^T^, are also reproduced by the last term in eq. (3.31) in 
our paper [46] (unfortunately additional factor 4 was missed there) 

An^^ ^^2 coth(ti?) 

8 R sinh^(ti?) 


Taking all this into account we arrive at the following low temperature behavior of the 
internal Casimir energy of a dilute dielectric ball 
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The relevant thermodynamic relations give the following low temperature expansions for 
free energy 
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( 5 . 42 ) 

( 5 . 43 ) 


The term in the free energy (5.42) does not give contribution to the Casimir force 
exerted on the surface of a dielectric ball, however, it proves to be important for insuring 
the positive entropy (5.43) at low temperatures. 

The range of applicability of the expansions (5.41), (5.42), and (5.43) can be roughly 
estimated in the following way. The curve S{T) dehned by eq. (5.43) monotonically goes 
up when the dimensionless temperature t = 2naT changes from 0 to t ~ 1.0. After that 
this curve sharply goes down to the negative values of S. It implies that eqs. (5.41) - (5.43) 
can be used in the region 0 < t < 1.0. The T^-term in eqs. (5.41) and (5.42) proves to be 
principal because it gives the hrst positive term in the low temperature expansion for the 
entropy (5.43). It is worth noting, that the exactly the same T^-term, but with opposite 
sign, arises in the high temperature asymptotics of free energy in the problem at hand (see 
eq. (4.30) in ref. [150]). 

For large temperature T we found [46] 

U(T)~^T, F(T) ~-^Tl\n(aT) - c], S(T) ~ ^l\n(aT) + c + 1], ( 5 . 44 ) 


where c is a constant [25,39] c = ln4 + 7 — 7/8. Analysis of eqs. (3.20) and (3.31) from 
the paper [46] shows that there are only exponentially suppressed corrections to the leading 
terms (5.44). 
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In the course of calculation of the thermodynamical functions of a dilute dielectric ball 
these functions were first obtained for a compact ball with the same velocity of light in¬ 
side and outside [46]. There we have derived the following exact, in the .^^-approximation, 
expression for the internal energy 
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In the small T region eq. (5.45) gives 
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Integration of the thermodynamic relation 


(7(r) = ^|/3F(T)], /3 = T-' 

enables one to get the free energy 

F(T) = -T j ^^dT + CT, 


(5.45) 


(5.46) 


(5.47) 


(5.48) 


where C is a constant. Substituting the asymptotics (5.46) into eq. (5.48) we obtain the 
respective free energy in the low temperature region 
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(5.49) 


Here the linear in T term CT has been dropped in view of the requirement that the entropy 
S{T) should vanish at T = 0 gives [36]. Indeed, the thermodynamic relation (7.14) gives in 
this limit 


^(0) = hm T-i {U{T) - F{T)) = C = 0 . 


(5.50) 


Hence, at low temperature the expansions both for the internal energy (5.46) and for the 
free energy (5.49) involve only even powers of the temperature beginning from T^. At zero 
temperature we have 


7/(0) = F(0) = Eball = 


5^^ 

327ra ’ 


(5.51) 


where i/baii is the exact, in the .^^-approximation, Casimir energy of a compact ball with the 
same velocity of light inside and outside the ball (see eq. (5.27)). 

The entropy in the problem at hand is obtained by differentiation of the free energy 
(5.49) (see eq. (7.14)) 
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The exact formula (5.45) leads to the following high temperature asymptotics of the 
internal energy 


U{T) ~ 



T 


oo. 


(5.53) 


Substituting this asymptotics into eq. (5.48) we arrive at the high temperature limit for 
the free energy F{T) 

F{T) ~ -^T [In(aT) + a], T ^ oo. (5.54) 

The constant a will be calculated in subsect. VIIC 2 (see eq. (7.42)) 


a = 7 + ln4 — 


5 


This constant has been found in ref. [36] too. The asymptotics (5.44) (5.53) and (5.54) 
contain terms the Planck constant, thus it is pure classical contributions (see also sect. VII). 

Summarizing we conclude that now there is a complete agreement between the results of 
calculation of the Casimir thermodynamic functions for a dilute dielectric ball carried out in 
the framework of two different approaches: by the mode summation method [42,46] and by 
perturbation theory for quantized electromagnetic held, when dielectric ball is considered as 
a perturbation in unbounded continuous surroundings [39]. 


VI. NON-SMOOTHNESS OF THE BOUNDARY AND DIVERGENCIES OF 

VACUUM ENERGY 

A. Semi-circular cylinder 

The spectrum of electromagnetic oscillations, as well as oscillations of other helds, is 
determined by the boundary form and by the conditions imposed on the held functions on 
the boundary. In view of this, one could anticipate that the divergencies in vacuum energy 
are also connected with the boundary geometry. In general it is true but this relation turns 
out to be very complicated and it is still far from being clear. We are going to show this 
by calculating the Casimir energy for a simple conhguration, namely, for a semi-circular 
cylindrical shall by making use of the zeta function technique. This shell is obtained by 
crossing an inhnite circular cylindrical shell by a plane passing through the symmetry axes of 
the cylinder. All the surfaces, including the inhnite cutting plane, are assumed to be perfectly 
conducting. Obviously it is sufficient to consider only a half of this conhguration (left or 
right) which we shall refer to as a semi-circular cylindrical shell or, for sake of shortening, as a 
semi-circular cylinder (see hg. 1). The internal boundary value problem for this conhguration 
is nothing else as a semi-cylindrical waveguide. In the theory of waveguides [119] it is well 
known that a semi-circular waveguide has the same eigenfrequencies as the cylindrical one 
but without degeneracy (without doubling) and safe for one frequency series (see below). 
Notwithstanding the very close spectra, the zeta function technique does not give a hnite 
result for a semi-circular cylinder unlike for a circular one (see subsect. IV C). 
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We start with considering the natural modes of electromagnetic held for circular and 
semi-circular cylinders. The construction of the solutions to the Maxwell equations with 
boundary conditions given on closed surfaces proves to be nontrivial problem. Mainly it 
is due to the vector character of the electromagnetic held [119,122,151]. In the case of 
cylindrical symmetry the electric E and magnetic H helds are expressed in terms of the 
electric (11') and magnetic (n") Hertz vectors having only one non-zero component 

n' = e,<l>(r,yp)e±'"^% (6.1) 

n" = e^T(r,(^)e^'^"". ( 6 . 2 ) 

Here the cylindrical coordinate system r, 99 , is used with axes directed along the cylinder 
axes. The common time-dependent factor is dropped. The scalar functions <h(r, 99 ) and 
'h(r, (p) are the eigenfunctions of the two-dimensional transverse Laplace operator and meet, 
respectively, the Dirichlet and Neumann conditions on the boundary ST 

(Vl+ 7'^)<I>(r,(^) = 0, $(r,(^)|gp = 0, (6.3) 
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( 6 . 6 ) 


First we consider a cylindrical shell. In this case the functions <h(r, ip) and T(r, ip) should 
be 27r-periodic in angular variable ip. As a result the Dirichlet boundary value problem (6.3) 
has the following unnormalized eigenfunctions (E-modes) 


^nm{r, ip) = (nip) 
cos 


( 1 ) 

nm 




, r < a, 
r), r > a, 


(6.7) 


where a is the cylinder radius, Jn(x) are the Bessel functions, II^^^(x) are the Hankel func¬ 
tions of the hrst kind, and 7 )^^, stand for the roots of the frequency equations 


Mjnma) = 0 , ( 6 - 8 ) 

n = 0,1, 2,..., m = 1, 2,... . 


For the Neumann boundary value problem (6.4) we have the Ff-modes 
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(6.9) 
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where 7 "^ and 7 "^ are the roots of the equations 

= 0 , ( 6 . 10 ) 

n = 0,1, 2,..., m = 1, 2,... . 

As usual, it is assumed that for r > a the eigenfunctions should satisfy the radiation condi¬ 
tion. 

It is important to note that each root 

7nm> 7nm> 7nm> 7nm> n > 1, m>l (6.11) 

is doubly degenerate since, according to eqs. (6.7), (6.9), because there are two eigenfunctions 
which are proportional to either sin(nv9) or cos{nif). The frequencies with n = 0 

7om> 7om) 7om) 7 om) ^ 1 ) 2 ,... ( 6 . 12 ) 

are independent on tp, and the degeneracy disappears. 

For given Hertz vectors 11' and 11" the electric and magnetic helds are constructed by 
the formulas 
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Jniin 


= 0 , 




'Wr 


E = VxVxn', H = —iu'V X n' (E-modes), 

E = iu'V X n", H = V X V X n" (i7-modes). (6.13) 


It has been proved [152] that the superposition of these modes gives the general solution 
to the Maxwell equations in the problem under consideration. An essential merit of using 
the Hertz polarization vectors is that in this approach the necessity to satisfy the gauge 
conditions does not arise. 

Now we turn to a waveguide which is obtained by cutting the inhnite cylindrical shell 
by a plane passing through the symmetry axes of the cylinder (see £g. 1). All the surfaces 
are assumed to be perfectly conducting. In this case the boundary value problems (6.3) and 
(6.4) for the Hertz electric (11') and magnetic (n") vectors have the following eigenfunctions 


^nm{r,ip) = sin(nv9) 

n = 1,2,... , 


r>a, 

m = 1, 2,... 


and 

'hnm(r, (p) = cos{np) 

n = 0,1,2,, m = 1,2,... . 


r<a, 

r>a. 


(6.14) 


(6.15) 


The frequencies 7 )^^, 7 ]^^) 7nm) 7 "^ are determined by the same equations (6.8) and 

(6.10) . However the new spectral problem has two essential distinctions: i) the frequencies 

(6.11) are now nondegenerate, and ii) two series of eigenfrequencies 

7L) m = l,2,... (6.16) 
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FIG. 1. The cross section of an infinite semi-circular cylindrical shell of radius a. All the surfaces 
(bold-faced lines) are assumed to be perfectly conducting. At the same time this picture presents 
the two-dimensional (plane) version of the problem under consideration, i.e., the semi-circular 
boundaries for massless fields defined on the plane. 


are absent. At first sight one could expect that such a change of the spectrum cannot 
influence drastically on the ultraviolet behavior of the relevant spectral density. However, 
as it will be shown below, the zeta function for a semi-circular cylinder, unlike for a circular 
one, does not provide a finite answer for the Casimir energy in the problem in question. 

In view of all above-mentioned the zeta function for electromagnetic field obeying the 
boundary conditions on the surface of the semi-circular cylinder is the sum of two zeta 
functions for scalar massless fields satisfying the Dirichlet and Neumann conditions on the 
lateral of this cylinder. 

By making use of the contour integration and analytical continuation discussed in sub¬ 
sect. IV C one can construct the spectral zeta functions for a semi-circular cylinder also. In 
ref. [153] it has been done first for scaler massless fields obeying, on internal and external 
surfaces of a semi-circle infinite cylinder, the Dirichlet and Neumann conditions. As a result 
the following expressions were obtained for relevant vacuum energies 
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(6.18) 


The sum of these formulae gives the Casimir energy of electromagnetic field for semi-circular 
cylinder 
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(6.19) 
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Pole singularities in eqs. (6.17), (6.18) and (6.19) imply that the zeta function technique does 
not provide a hnite value of vacuum energy for this conhguration and further renormalization 
is needed. 

The same situation takes place also in the two-dimensional version of the problem under 
study, i.e., when a semi-circle is considered on a plane. In fact, the relevant zeta functions 
can be obtained in a straightforward way by making use the relation (4.76). It gives [153] 
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( 6 . 20 ) 


Ef-cir = (-0 = ^ (-0.237103 - 0.0062 ^ J . (6.21) 

The sum of these energies is also inhnite. 

Keeping in mind that the zeta function regularization does nor supply a hnite value for 
the Casimir energy also in the case of spheres in even dimensional spaces (see subsect. IV C) 
we have thus a quite long list of boundary conditions for which the vacuum energy cannot 
be calculated by making use of the known mathematical methods. It is worth elucidating 
the geometrical peculiarities of the boundaries that are responsible for this failure. In order 
to do this we have to address the corresponding heat kernel coefficients. 


B. Mutual cancellation of divergencies in vacuum energy 

The relevant analysis of the heat kernel coefficients (up to B 5 / 2 ) has been accomplished 
in ref. [116]. We are interested in the coefficients B 2 and i? 3 / 2 - As it has been noted 
in subsect. IIIC, the nonzero value of B 2 {B 3 / 2 ) implies that it is impossible to derive a 
hnite value for the vacuum energy in the problem at hand, when the dimension of the 
conhguration space d = 3 {d = 2). From the asymptotic expansion (3.25) it follows that 
the heat kernel coefficients introduced in this way are the same for the three-dimensional 
cylindrical-like boundaries and for the corresponding plane problem obtained by crossing the 
former boundary by a transverse plane. Certainly, it is correct only for the Laplace operator 
—A. In fact, the respective eigenvalues in these two problems are related by 

Xn{d = 3) = k'^ + Xn{d = 2), 0 < k < 00 . (6.22) 


Hence 

/ OO JU 1 

e-^^^K,= 2 {t) = -=iF,=2(t). (6.23) 

-00 27r V47rt 

Taking into account the dehnition (3.25) one can easily deduce form (6.23) that the heat 
kernel coefficients Bn /2 in the two eigenvalue boundary value problems mentioned above, 
are equal. 
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TABLE I. The contribution of different parts of boundary to heat kernel coefficients; D and 
N stands for the Dirichlet and Neumann boundary conditions; the upper (lower) sign refers to 
internal (external) region. 
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The heat kernel coefficients we are interested in are contributed by the following pecu¬ 
liarities of the boundary (see £g. 1 ): the curvature of semi-circle 1-2 and four right-angled 
corners at points 1 and 2. It is worthy to remind that we are considering both internal (-Din) 
and external (Dgx) regions. All these contributions are presented in table I. The contribu¬ 
tions of the curvature of the semi-circle 1-2 to the coefficients B 2 for internal and external 
regions are of the same absolute value but they have opposite signs. As a result they are 
mutually cancelled in the net B 2 coefficient. Unlike this the contributions to B 2 due to four 
right-angled corners at the points 1 and 2 are added 
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(6.24) 


This fact has been noted at first time in ref. [154]. In the case of the coefficient D 3/2 the 
situation is opposite, i.e., the contributions of the corners from internal and external regions 
are mutually cancelled, while the contributions of the curvature of the arc 1-2 from internal 
and external regions are added 


,D _ rN _ p STTy/j _ STTy/j _ STTy/j 

64D 32D’ 64D 32R ’ 2/2 • 


(6.25) 


Different geometrical origins of the zeta function failure to provide a finite value of the 
vacuum energy in the two- and three-dimensional versions of the boundary value problem 
in question probably imply the impossibility of obtaining a finite and unique value of this 
quantity by taking advantage of the atomic structure of the boundary [155] or its quantum 
fluctuations [156]. It is clear because any physical reason of the Casimir energy divergences 
should be valid simultaneously in the two- and three-dimensional versions of the boundary 
configuration under consideration. 
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The alternating sign changes of contributions due to curvature and corners to heat kernel 
coefficients are presumably true for all the coefficients with n=l, 3/2, 2,..., at least 
it is the case for i?i, -B 3 / 2 , -B 2 and -B 5/2 (see table I and also ref. [116]). The analogous 
situation takes place for spheres in spaces of even and odd dimensions, namely, the curvature 
contributions to respective heat kernel coefficients are mutually cancelled for odd dimension 
of ambient space and they are added in spaces of even dimension (see subsect. IV C). In 
the former case the zeta regularization gives a hnite value for vacuum energy, in the latter 
situation it is not the case, and further renormalization of vacuum energy is needed in the 
problem at hand. 

The corners at points 1 and 2 of the boundary considered here (see £g. 1) are obtained 
by intersection of a straight line and an arc of a circle possessing nonzero curvature ( 1 /a). 
They contribute to all the heat kernel coefficients starting with n = 1. These corners should 
be distinguished from those formed by crossing two straight lines, for example, from the 
corners of a rectangle. Such corners contribute only to the coefficient Bi, the contribution 
being the same for Dirichlet and Neumann boundary conditions 

c(a) = (6.26) 


where a is the angle of the corner. Such corner singularity of the boundary proves to be 
very close to the conical singularity, which is important in many areas of mathematical 
physics [157,158]. Lately it has been investigated in connection with studies of quantum 
fields on the background of black holes [159] and cosmic strings [160,161]. 

This general assertion concerning the corner contribution to the heat kernel coefficients 
can be illustrated by a known heat kernel expansion for a rectangle with sides a and b 


m=l n=l 



ab 

Ant 


a + b 
A\/ni 


1 

4 


ES, 





+ ES 


(6.27) 


where ES denotes the exponentially small corrections as t ^ +0 (see, for example, [162]). 
Here the scalar operator —A with Dirichlet boundary conditions is considered. For a rect¬ 
angle the coefficient Bi is obviously equal to the contributions of four right-angled corners. 
Indeed, the third term in the expansion (6.27) can be represented in the form 


1 

4 




V2), 


where c{a) is given in eq. (6.26). Besides these corners the boundary of a rectangle has no 
other singularities, therefore the heat kernel coefficients Bn with n = 3/2, 2, ... vanish in 
this problem. 

These rules for obtaining the heat kernel coefficients are directly generalized to an arbi¬ 
trary polygon with the angles ctj. The first two coefficients Bq and H 1/2 are defined by eq. 
(3.26), where V is the area of the polygon and S is its perimeter. The third coefficient Bi 
is equal to the sum of the contributions due to the angles aj 
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(6.28) 


The rest of the coefficients 3^-, n > 1 vanish both for internal and external regions. In 
particular, it implies that the zeta function technique should provide a finite value of the 
Casimir energy for a polygon on a plane (i? 3/2 = 0) and for a cylindrical generalization of 
the polygon spectral problem (i ?2 = 0). These subjects have been discussed earlier in pa¬ 
pers [163]. The calculations of vacuum energy have been accomplished only for parallelepiped 
geometries [164-173] and with allowance for only internal region. For such conhgurations 
the spectrum of the Laplace operator is known exactly and the zeta regularization works 
well. For example, the Casimir energy for the cube with sides a is negative 

ri ri 1 ^ 

Ecnhe = ^ \/v? +w? + kP‘ ~ — ^- . (6.29) 

2a a 

n,m,fe=l 

However Milton [59] comments these results in the following way: here ’’divergences occur 
which cannot be legitimately removed, which nevertheless are artihcially removed by zeta 
function methods. It is the view of the author that such finite results are without meaning”. 
With regard to external region the wave equation is not separable outside a cube or a 
rectangular solid. 

Recently the image method has been extended to a set of geometries with planar bound¬ 
aries and without rectangular angles [174]. In this approach the Casimir energy of a mass¬ 
less scalar held was calculated inside a triangle with angles proportional to 7r/iV, N > 3. 
However the dual expression for the vacuum energy is not ready for numerical estimation. 
Further this technique was applied for calculating the Casimir energy of a conical wedge and 
a conical cavity [175]. 


VII. HIGH TEMPERATURE BEHAVIOR OF THE VACUUM ENERGY 

The inhuence of temperature on the Casimir effect was an important topic since its first 
experimental demonstration [2] which had been done at room temperature. It was first 
shown in ref. [176] that the temperature inhuence was just below what had been measured 
(see, also, paper [177]). Now it is expected that the temperature contributions will be seen 
in the upcoming series of experiments [10]. The temperature dependence of the vacuum 
energy proves to be important practically in all problems where the Casimir ehect is taken 
into account, for example, in hadron physics (quark deconhnement in the framework of bag 
models [178]) and in cosmology [179]. 

The Casimir calculations at finite temperature is a nontrivial problem specihcally for 
boundary conditions with nonzero curvature. Investigation of the high temperature limit, 
i.e., , the classical limit, in this problem is of independent interest [180]. For this goal a 
powerful method of the zeta function technique and the heat kernel expansion can be used. 
It is important that for obtaining the high temperature asymptotics of the thermodynamic 
characteristics it is sufficient to know the heat kernel coefficients and the determinant for 
the spatial part of the operator governing the held dynamics. This is an essential merit of 
this approach. 
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A. Heat kernel coefficients and high temperatnre expansions 


In quantum field theory, finite temperature effects can be described at equilibrium in 
the Matsubara formalism by imposing periodic (resp. antiperiodic for fermions) boundary 
conditions in the imaginary time coordinate [181]. Let the dynamics of quantum field (p(t, x) 
be defined, as before, by eq. (3.3). The Helmholtz free energy F of this field is determined 
by the functional integral [181,182] 


,-l3F _ 


= / Dip exp 


dt dx ip*(t, x)LTip{t, x] 


(7.1) 


Here (3 is inverse temperature f3 = T~^. For simplicity the Boltzmann constant k-Q is 
assumed to be equal to 1. Therefore the temperature T is measured in energy units. For the 
Euclidean time and for the Euclidean version of the field ip{t, x) we use the former notations. 
The operator L in eq. (7.1) is the Euclidean version of the full differential operator in the 
field equation (3.3) 


with the eigenfunctions 


L'j' — L — 


dD 




(7.2) 


(7.3) 


and the eigenvalues 

k = {m,n}, = + (7.4) 

where Qm = 2nmT/h, m = 0,±1,±2,... are the Matsubara frequencies, and ojn are 
defined in eq. (3.4). 

Functional integration in eq. (7.1) yields 

= (dethT)-^/^. (7.5) 

The determinant of the operator is expressed in terms of its spectral zeta function 

OO 

C,t(s) = = Y. + ( 7 . 6 ) 

k m=—oo n 

Indeed 

det L(r = A^, ln(det Lt) = ^ In Aj. (7.7) 

k k 

Differentiation of the zeta function Ct{s) at the point s = 0 gives 

Ct(0) = -E'“T (7.8) 

k 
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/,From eqs. (7.6), (7.7) and (7.8) we deduce finally 

f = -|cT(0). (7.9) 

The characteristics of the quantum field (3.3) at zero temperature are determined by the 
zeta function (C(s) associated with the operator L (see eq. (3.9)). From the mathematical 
point of view the zeta function (C(s) corresponding to the space part of the full wave equation 
(3.3) is, undoubtedly, a simpler object than the complete zeta function Cr('S) because the 
definition (7.6) involves an additional sum over the Matsubara frequencies. Here a natural 
question arises whether one can gain knowledge of the quantum field at nonzero tempera¬ 
ture possessing only the zeta function C('S)- In ref. [183] it was shown that proceeding from 
the zeta function (C(s) one can deduce the high temperature asymptotics of the thermody¬ 
namic functions such as Helmholtz free energy, internal energy, and entropy. Let us remind 
briefly the derivation of these asymptotics. By making use of the formula (3.24) the zeta 
function (3.23) can be represented in the form 


Ct{s) 



(7.10) 


The term with m = 0 in this formula gives the zeta function (3.9). In the remaining terms 
we substitute the heat kernel K{t) of the operator L by its asymptotic expansion (3.25). As 
a result we arrive at the following asymptotic representation for the complete zeta function 
Ct(-s) 


Ct(s) ^ C(s) 


2 

(27r)3/2 


n=0,l/2,... 


h y"^+^'^F(s-3/2 + n) 

f(i) 


Cr(2s -h 2n - 3), 


(7.11) 


where Cr('S) is the Riemann zeta function (3.11). Taking the derivative of the right hand 
side of eq. (7.11) at the point s = 0 and substituting the result into eq. (7.9) one obtains 
the high temperature expansion for the free energy 


T TT^ 

F(T).--C'(0)+il.^- 
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24 h ^ (47r)3/2 
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Idvr^ 

-rE 
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(47r)3/2 \27iT 
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2n—3 


5/2 




(47r)3/2 24T 


F(n — 3/2) CR(2ri — 3), T —> cx). 


(7.12) 


Here 7 is the Euler constant. The argument of the logarithms in expansion (7.12) are 
dimensional, but upon collecting similar terms with account for the logarithmic ones in 
C'( 0 ) it is easy to see that finally the logarithm function has a dimensionless argument, at 
least for Bo = 0. It is worth noting that the auxiliary evolution variable r in eq. (3.25) has 
the dimension [length]^. 
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The asymptotic expansions for the internal energy U{T) and the entropy S{T) are de- 
dnced from eq. (7.12) employing the thermodynamic relations 


d 


U{T) = -T^— (T-'F(T)) , 
S(T) = T-^ {U{T) - F{T)) = 


dF 

’&r' 


(7.13) 

(7.14) 


Snbstituting the expansion (7.12) into eqs. (7.13) and (7.14) one arrives at the asymptotics 


U{T) ~ Bo 
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(47r)3/2 12 T 
r(n-l/2)CR(2n-3), 
T 


(7.15) 
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/ k \ 2n—3 

E ^» ( ^ ) (« - 2) r(n - 3/2) CR(2n - 3). (7.16) 

n>3 ^ 2^ 


In eq. (7.15) the term proportional to B 2 contains the logarithm of dimensional qnantity: 
[h/T] = [time]“^. This is the result of the arbitrariness arising in from the ultraviolet 
divergences in the case of i ?2 4 0 (see ref. [146] for a more detailed discussion). Unlike this 
situation, collecting the logarithm functions in the i? 3 / 2 -term and in C^(0) in eq. (7.16) leads 
to a dimensionless argument of the logarithm in the hnal expression. 

It is worth noting that the zeta determinant of the operator L, i.e., C^(0)) does not enter 
the asymptotic expansion for the internal energy (7.15). Therefore this high temperature 
expansion is completely dehned only by the heat kernel coefficients. In view of this, the 
hrst term in the asymptotics of the free energy in eq. (7.12) is referred to as a pure entropic 
contribution. Its physical origin is till now not elucidated. 


B. Perfectly conducting parallel plates in vacuum 

First we demonstrate the application of the high temperature expansions (7.12), (7.15) 
and (7.16) to a simple problem of electromagnetic held conhned between two perfectly 
conducting parallel plates in vacuum. 

Substituting eq. (4.4) into eq. (3.27) we obtain for perfectly conducting parallel plates 
only one nonzero coefficient Bq 


Bo = 



(7,17) 
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where V = L^Ly a is the volume of the space bounded by the plates.^ In subsect. IV A it 
was noted that for parallel conducting plates electromagnetic field is reduced to two massless 
scalar fields obeying Dirichlet and Neumann conditions on internal surface of the plates. As 
a result we have the multiplier 2 in eq. (7.17) and vanishing coefficient B 1/2 (see eq. (3.26)). 
It should be noted here that we are considering only electromagnetic held conhned between 
the plates and do not take into account that held outside the plates. 

^From eqs. (7.15) and (7.17) it follows that the density of internal energy has the following 
high temperature asymptotics 


U(T) 

V 




4-T^ 


T 


CX), 


(7.18) 


where a is the Stefan-Boltzmann constant 


60 c2 ■ 


(7.19) 


Recall that in our formulae we put /cb = 1, that is, the temperature is measured in energy 
units. The transition to degrees is performed by the substitution T k^T. 

When calculating the high temperature asymptotics of the free energy (7.12) and the 
entropy (7.16) one needs to derive C^(0) fhe zeta function (4.4). Keeping in mind that 
(Cb(— 2) = 0 it is convenient to use here the Riemann rehection formula [112] 


2^ "r(s)CR(s) cos(7rs/2) = 7r^CR(l - s) 


(7.20) 


which yields 

^From here we deduce 


(7.21) 


C(0) = ^Ck(3) = j^Ck(3). 


(7.22) 


Insertion of eqs. (7.17) and (7.22) into eq. (7.12) gives the following high temperature be¬ 
havior for the density of free energy 


V c3h3 90‘ 

As was noted above, we are considering only electromagnetic held between the plates. There¬ 
fore when calculating the Casimir forces one should drop the last term in eq. (7.23) since its 
contribution is cancelled by the pressure of the black body radiation on the outward surfaces 



^For obtaining the vanishing ili /2 coefficient it is important to take into account the second term 
in eq. (4.4) which depends on the photon mass 
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of the plates. As a result the high temperature asymptotics of the Casimir force, per unit 
surface area, attracting two perfectly conducting plates in vacuum is 


T ^ 


T 

4:7ia^ 


Cr(3). 


(7.24) 


Usually in the Casimir calculations the contribution of the free black body radiation is 
subtracted from the very beginning [67]. 

It is interesting to note that the Casimir force (7.24) and the hrst term on the right 
hand side of eq. (7.23) are pure classical quantities because they do not involve the Planck 
constant h. These classical asymptotics seem to be derivable without appealing to the notion 
of quantized electromagnetic held. The classical limit of the theory of the Casimir effect is 
discussed in a recent paper [180]. 

Employing eqs. (7.14) and (7.23) one arrives at the high temperature behavior of the 
entropy density 


S{T) ^Cr(3) 2T^7t^ 

V 87ra3 ^ 45c3h3‘ 


(7.25) 


Corrections to eqs. (7.18), (7.23) and (7.25) are exponentially small. 


C. Sphere 

We consider electromagnetic held snbjected to three types of bonndary conditions on 
the snrface of a sphere: i) an inhnitely thin and perfectly condncting spherical shell; ii) 
the surface of a sphere delimits two material media with the same velocity of light; iii) a 
dielectric ball placed in unbonnded dielectric medium. In order to obtain the heat kernel 
coefficients determining the high temperatnre asymptotics (7.12), (7.15) and (7.16) it is 
convenient to use the explicit representation of the relevant spectral zeta fnnctions in terms 
of the Riemann zeta function. These formulae were derived in subsect. IV B by taking into 
account the hrst two terms of the uniform asymptotic expansion for the product of the 
modihed Bessel functions Iy{iJz) Ky{iJz). 


1. Perfectly conducting spherical shell 

The corresponding spectral zeta function is given in eq. (4.33). The terms omitted in 
this eqnation are of the form 

gfc(s) [(22(^+")+^-1) Cr(2A; + 2s + 1)-22(^+^)+i] , A; = 2, 3,4,..., (7.26) 

where qk{s) stand for some polynomials in s. 

Analysis of eqs. (4.33) and (4.32) shows that the zeta fnnction for a perfectly condncting 
spherical shell enables one to hnd the exact values of the hrst six heat kernel coefficients, 
namely: 

Ro = 0, Ri/ 2 = 0, Ri = 0, R3/2 = 27r^/^, R 2 = 0, 55/2 = ——(7.27) 

2U 
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Taking into account the structure of the omitted terms (7.26) it is easy to see that 


= 0, j=3,4,5,... . 


(7.28) 


Having obtained the heat kernel coefficients (7.27) and (7.28) we are in position to construct 
the high temperature asymptotics of the internal energy of electromagnetic held by making 
use of eq. (7.15) 


U{T) 




1 

1920 T 


+ C>(T-3). 


(7.29) 


Applying the technique developed in ref. [184] more terms to this expansion can be easily 
added. 

In order to write the asymptotic expansions (7.12) and (7.16) the derivative of the zeta 
function at the point s = 0 should be calculated. Equation (4.33) gives an approximate 
value for C^O) 

Csheii(O) ^ 2 16 ~ 8 2 c ^ 0.38265 + - In - . (7.30) 

The terms omitted in (4.33) will render precise only the hrst term in the hnal form of this 
expression, while the second term (1/2) ln(a/c) will not change. The exact value of Cshell(O) 
is calculated in Appendix B 1 

Csheii(O) = ^- | + ^ln2 + 6 Cr(- 1) + (^-^ + ^ In ^ + In 2 + 

= 0.38429 +-In-. (7.31) 

2 c 

It is worth noting that the expression in the round parentheses, being multiplied by is 
exactly the value of Cbaii(O) ^ compact ball with continuous velocity of light on its surface 
(see eq. (7.41) in the next subsect.). As a result we have the following high temperature 
asymptotics of the free energy and the entropy in the problem in question 

F(T)~-- fin ^+ 0.76858^ - f — 1 ^ — + 0{T-^), (7.32) 

^ ^ A \ he ) \a ) 3840T ^ ^ ^ 

SiT) ~ 0.44215 + - In ^ ^ f—] + 0(T -^). (7.33) 

^ ^ 4 he 3840 [aTJ ^ ^ ^ ^ 


The expression (7.32) exactly reproduces the asymptotics obtained in ref. [15] (1978) by 
making use of the multiple scattering technique (see eq. (8.39) in that paper). We have not 
calculated the coefficient i? 7 / 2 , therefore we do not know the sign of the T“^-correction in 
(7.32). In ref. [15] it is noted that this term is negative. 

In eqs. (7.29), (7.32) and (7.33) the large expansion parameter is actually a dimensionless 
‘temperature’ r = aT/{hc). Therefore the same formulae describe the behavior of the 
thermodynamic functions when a ^ oo and temperature T is fixed. 

The high temperature asymptotics of the thermodynamic functions derived by making 
use of the general expansions (7.12), (7.15) and (7.16) contain the terms independent of the 
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Planck constant h or, in other words, classical contributions (see eqs. (7.29), (7.32) and 
(7.33)). This is also true for the high temperature limit of the Casimir force calculated per 
unit area of a sphere 


J^(T) ~ 


1 dF{T) 
47ri?2 da 


T 

IhTra^ 


hcV 1 1 

a ) 47ia^ 1920 T 


+ C>(T-3). 


(7.34) 


The leading classical term in the asymptotics (7.34) describes the Casimir force that 
seeks to expand the sphere. The quantum correction in this formula stands for the Casimir 
pressure exerted on the sphere surface. 

In eqs. (7.29), (7.32) and (7.33) the Stefan-Boltzmann terms proportional to are absent 
because the contribution of the Minkowski space was subtracted from the very beginning in 
our calculations (see eq. (4.12)). As a result we obtain the vanishing heat kernel coefficient 
Bq which, in general case, is determined by the volume of the system under study (see eq. 
(7.17)). Therefore our results describe only the deviation from the Stefan-Boltzmann law 
caused by the perfectly conducting sphere. 

The vanishing of the coefficients B 1/2 and Bi in the problem at hand can be explained 
by taking into account the general properties of the heat kernel coefficients [60] and by 
making use of the results obtained in ref. [184]. As known [122] the solutions to the Maxwell 
equations with allowance for a perfectly conducting sphere are expressed in terms of the two 
scalar functions that satisfy the Laplace equation with the Dirichlet and Robin boundary 
conditions on internal and external surfaces of the sphere. In view of this one can write 

+ B^_ + B^_^_ + B^_, n = 1/2, 1, ... , (7.35) 


where the subscript plus (minus) corresponds to internal (external) region and the rest 
notations are obvious. In ref. [184] it was found that 

51 / 2 + = -27r3/2a2 = ^V2+ = = B^/ 2 -^ 

= ±^a, Bf^ = (7.36) 

As a result we have for electromagnetic held inside and outside a perfectly conducting sphere 


Ri/2 = = 0 . (7.37) 

Having calculated the corrections to the Stefan-Boltzmann law one should naturally 
discuss the possibility of their detection. The ratio of the leading term in eq. (7.29) to the 
internal energy of black body radiation in unbounded space given by the Stefan-Boltzmann 
law (7.18) is proportional to r“^. Already for r ~ 10 the corrections prove to be of order 
10“^. The same value of r can be reached by varying the scale of length a in the problem 
under consideration or by respective choice of the temperature T. Keeping in mind the 
value of the conversion coefficient ch = 197.326MeV ■ fm = 0.229K ■ cm [185] we obtain the 
following estimations. For a ~ 10“^^cm (a typical hadron size) the temperature T should 
satisfy the inequality T 3> 200MeV in order to apply the asymptotics found. For a ~ 1cm 
we have T ^ 0.229K and for a ~ 7 ■ lO^^cm (radius of the Sun) the range of applicability 
of the asymptotics at hand extends practically to any temperature value T ^ 10“^°K. Here 
we shall not go into the details of a concrete experimental equipment that enables one to 
observe the calculated corrections to the Stefan-Boltzmann law conhning ourselves to the 
estimations presented above. 
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2. Compact hall with equal velocities of light inside and outside 


The spectral zeta function for this conhguration is given in eqs. (4.27) and (4.25). It 
affords the exact heat kernel coefficients up to i? 5/2 

Bq = 0, i?i/2 = 0, Bi = 0, i?3/2 = i?2 = 0, 

^5/2 = ^^%7rv/ip(-l) = 0. (7.38) 

CL 

With allowance of the structure of the omitted terms in eq. (4.27) we can again deduce that 

B, = 0, j = 3,4,5,.... (7.39) 

Substitution of these coefficients into eq. (7.15) gives the following high temperature behavior 
of the internal energy in the problem under consideration 

U{T) ^ ^ + 0(T-3). (7.40) 

The value of Cbaii(O) calculated in the Appendix B2 

C'(0) = ^ 7 + '“ 2 + l) = (o.35676 + 1 In ^ j . (7.41) 

It is this value that is supplied by eq. (4.24) with allowance for that p(—1) = 0. 

By making use of eqs. (7.12), (7.38) and (7.41) we deduce the high temperature asymp¬ 
totics for free energy 

F(r) = -yl (^7 + i„4-5j +^rinl^ + o(r-3), 

= -£2 ^0.71352 + ^ T In — + 0(T-^). (7.42) 

4 4 aT 

The entropy in the present case has the following high temperature behavior 

s(T) = ^ (^i + 7 + i„4-5-i„l£^+o(r-‘), 

= 4 (1.71352 - In 4^) + 0(T-‘). (7.43) 

4 \ aT J 

The asymptotics (7.40) and (7.42) completely coincide with the analogous formulae obtained 
in subsect. VB by the mode summation method combined with the addition theorem for 
the Bessel functions (see also ref. [36]). 

The exact expression for the internal energy in the problem at hand (see eq. (5.45)) gives 
only exponentially suppressed corrections to the leading term (7.40) 

U{T) ~ [1 + 2{4£^ + 4t + l)e-^*] . (7.44) 
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The asymptotics (7.44) implies in particular that in reality in eq. (7.40) there are no cor¬ 
rections proportional to the inverse powers of the temperature T. From here it follows 
immediately that all the heat kernel coefficients with integer and half integer numbers equal 
or greater than 3 should vanish 

B, = 0, j= 3, 7/2, 4, 5/2, 5, ... (7.45) 

(compare with eq. (7.39)). In view of this the sign denoting the omitted terms in eqs. (7.40), 
(7.42) and (7.43) should be substituted by (9(e“®™'^). 


3. Dielectric ball in unbounded dielectric medium 


The zeta function for electromagnetic held in the background of a pure dielectric ball 
(/ii = /i 2 = 1, 7 ^ £ 2 ) has not been obtained in an explicit form. In ref. [146] the heat 

kernel coefficients up to B 2 in this problem were found. Here we use the results of this paper 
conhning ourselves to the An^-approximation, where An = ni — n 2 = nin 2 (c 2 — ci)/c ~ 
(c 2 —Cl)/c, Ui and Cj are the refractive index and the velocity of light inside (i = 1 ) and outside 
{i = 2 ) the ball, and c is the velocity of light in the vacuum: n^ = Ci = c/nj, i = 1 , 2 . 

It is assumed that Ci and C 2 differ from c slightly, therefore C 2 —Ci and An are small quantities. 
In view of this we have 

Ho = ~ (An -I- 2 An^), 

Hi/ 2 = -27r3Aa^ ^ _4^3/2 ^^^ 2 ^ 

Hi ~ 0, H 3/2 = ^ TT^/^An^, H 2 ~ 0. (7.46) 

The coefficients Hi and H 2 equal zero only in the An^-approximation considered here. In 
the general case they contain terms proportional to An^, where k > 3. 

Allowance for one more term in the uniform asymptotic expansion of the modihed Bessel 
functions, as compared with the calculations in ref. [146], gives the next heat kernel coefficient 


(47r)3/2 2688 


(7.47) 


Correcting the mistake made in [186] we state that this coefficient has no contributions 
proportional to An^, and in the An^-approximation one has to put 


H 5/2 ^ 0 . (7.48) 

Making use of the technique developed in ref. [187,188] one obtains the following expres¬ 
sion for the derivative of the zeta function for a pure dielectric ball at the point s = 0 (see 
ref. [150]) 

C'(0) = ^(-g + ln j + ln4 + i), (7,49) 
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Before turning to the construction of the high temperature asymptotics in the problem 
at hand by making use of the general formulae (7.12), (7.15) and (7.16) a physical remark 
should be done. When considering the electromagnetic held in the background of a dielectric 
body in the formalism of quantum electrodynamics of continuous media, as a matter of 
fact one is dealing with a system consisting of two objects: electromagnetic held plus a 
continuous dielectric body. It is important that this body is described (phenomenologically) 
only by respective permittivity without introducing into the Hamiltonian special additional 
dynamical variables. As a result the zeta function and the relevant heat kernel coefficients 
calculated in this formalism also describe both electromagnetic held and dielectric body. 
When we are interested in the Casimir thermodynamic functions in such problems we have 
obviously to separate in the general expressions the contributions due to the dielectric body 
itself [189]. 

Let us turn to such separation procedure in the high temperature asymptotics for a 
dielectric ball. Following the reasoning of refs. [47,148] we divide the Helmholtz free energy 
of a material body with volume V and the surface area S into the parts 

F = Vf + Sa + Fc^, (7.50) 

where / is the free energy of a unit volume of a ball, a denotes the surface tension, and 
Fcas is refereed to as the Casimir free energy of electromagnetic held connected with this 
body and having the temperature T. In this way we obtain the following high temperature 
behavior of the free energy F(T) in the problem at hand 

where Bq and B 1/2 are dehned in eq. (7.46) and 

FcUT) ~ T (in L4 + 0 + o(T-^). (7.52) 

The high temperature asymptotics for the Casimir internal energy and for the Casimir 
entropy can be derived by making use of the respective thermodynamical relations (7.15), 

(7,16) 

77cas(T)^^T + (!l(T-2), (7.53) 

Sc.,(T) c. ^ (1 + 7 + In + o(T-^). (7.54) 

It is worth comparing these results with analogous asymptotics obtained by diherent 
methods. In our paper [46] at the beginning of calculations the hrst term of expansion 
of internal energy (7.53) was derived. The subsequent integration of the thermodynamic 
relation (7.13) gave the correct coefficient of the logarithmic term in the asymptotics of free 
energy (7.52) (see eq. (5.44)). In paper [47] Barton managed to deduce the asymptotics 
(7.52) - (7.54). One should keep in mind that our parameter An corresponds to 27ran in 
the notations of ref. [39]. 
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The asymptotics (7.52)-(7.54) contain the a-independent terms. As far as we know the 
physical meaning of such terms remains unclear. 

Preliminary analysis of a complete expression for the internal energy of a dielectric ball 
(see eqs. (3.20) and (3.31) in ref. [46]) shows that probably there are only exponentially 
suppressed corrections to the leading term (7.53). In that case in addition to eq. (7.48) all the 
heat kernel coefficients with number greater than 3 should vanish in the An^-approximation. 


D. Cylinder 

The calculation of the vacuum energy of electromagnetic held with boundary conditions 
dehned on a cylinder, to say nothing of the temperature corrections, turned out to be 
technically a more involved problem than the analogous one for a sphere. Therefore the 
Casimir problem for a cylinder has been considered only in a few papers [15,19-21,26,49,124]. 
We again examine three cases: i) perfectly conducting cylindrical shell; ii) solid cylinder with 
Cl = C 2 ; iii) dielectric cylinder when ci 7 ^ C 2 . Here we shall use the results of our previous 
papers [49,124]. 


1. Perfectly conducting cylindrical shell 


The zeta function for this problem is given by eqs. (4.70), (4.55), (4.61) and (4.71). 

The function Zi{s) is dehned in the strip —3/2 < Res < 1/2, while the functions 
^ 2 ( 5 ) and ^’ 3 ( 5 ) are analytic functions in the whole complex plane s except for the points, 
where r(s) and Cr('S) have simple poles. In order to hnd the heat kernel coefficients Bq, 
Ri/ 2 , and Bi trough the relation (3.27) one needs the zeta function dehned in the region 
1/2 + e < Res < 3/2 + e with e being a positive inhnitesimal. However in this region 
eq. (4.55) is not applicable directly due to the bad behavior of the integral at the upper 
limit. In the most simple way we can overcome this difficulty as in the case of perfectly 
conducting plates by introducing the photon mass p at the very beginning of the calculation 
and making then the analytic continuation of the zeta function to the points s = 1/2, 1, 3/2. 
Upon taking the residua at these points one should put p = 0. 

With regard to all this and using the relation (3.27) we hnd the heat kernel coefficients 


Bq — 0 , Bi/2 — 0 , Bi — 0 , B2 — 0 . 


(7.55) 


The vanishing heat kernel coefficient B 2 implies that the zeta regularization gives a hnite 
value for the vacuum energy in the problem at hand (see eq. (4.72) and ref. [19]). The 
coefficient B 3/2 is determined by the function Z 2 {s) only (see eq. (4.61)) 


B3/2 _ 3 

( 471 -) 3/2 54 Q 


(7.56) 


The coefficient R 5/2 is dehned by the function Z^i^s) given in eq. (4.71) 

-B5/2 _ 153 
(47r)3/2 “ 81^^‘ 


(7.57) 



The calculation of the next heat kernel coefficients -B 3 , -B 7 / 2 ,... would demand a knowledge 
of the additional terms in the expansion of the spectral zeta function in the problem under 
consideration in terms of the Riemann zeta function. These terms are proportional to 
Cr(2A; + 2s + 1) with k = 2,3,..., and may be obtained employing the technique presented 
in subsect. IV C. Analyzing the position of poles for these Riemann zeta functions it is easy 
to show that, as well as in the spherical case, we have 


R, =0, j = 3, 4, 5.... 


The zeta determinant entering the high temperature asymptotics of Hlemholtz free en¬ 
ergy (7.12) and entropy (7.16) is calculated in the Appendix Cl 


C'(o) 


0.45847 3 , a 

-^-In—. 

a 32 a 2 c 


(7.58) 


Now we are able to construct the high temperature expansions of the thermodynamic 
functions in the problem under consideration. For the free energy we have 


T 

F{T) ~ -0.22924- 
a 


3T , aT 

- In —— 

64 a 2 he 


51 

65536 


+ 0{T-^). 


(7.59) 


When comparing eq. (7.59) with results of other authors one should remember that all the 
thermodynamic quantities that we obtained in this section are related to a cylinder of unit 
length. The high temperature asymptotics of the electromagnetic free energy in presence of 
perfectly conducting cylindrical shell was investigated in ref. [15]. To make the comparison 
handy we rewrite their result as follows 


rri Q rri rri 

F(T) ~ -0.10362-—In (7.60) 

^ ^ a 64a 2hc ^ ^ 

The discrepancy between the terms linear in T in eqs. (7.59) and (7.60) is due to the double 
scattering approximation used in ref. [15] (see also the next subsection). Our approach 
provides an opportunity to calculate the exact value of this term (see eq. (7.59)). 

And hnally, making use of the general formulae (7.15) and (7.16) we derive 


U{T) ~ 

R(T)^ 


3T 153 c^h^ 
^ ~ 98304 

0.27612 3 , 

-^-In 

a 64 a 


■ + 0(T-3), 

RT 153 
2hc ~ 196608 a3T2 


+ 0{T-^). 


(7.61) 

(7.62) 


2. Compact cylinder with ci = C 2 and with ci / C 2 

Here we consider the boundary conditions for electromagnetic held of two types: i) a 
compact inhnite cylinder with uniform velocity of light on its lateral surface, ii) a pure 
dielectric cylinder with ci 7 ^ C 2 . 

The zeta function for the former conhguration is given by eqs. (4.60), (4.56), (4.58) and 
(4.59). It enables one to hnd the following heat kernel coefficients 
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Bq — 0, Bi /2 — 0, Bi — 0, -B3/2 — ; B2 — 0, 

^5/2 = 5, = 0, j = 3,4,5,... . (7.63) 

(Jj 

As before we are considering the ^^-approximation. The heat kernel coefficients (7.63) lead 
to the following high temperature behavior of the internal energy in the problem at hand 


u{T) 


64 i? V* 512i?2T2j 

The corresponding zeta determinant is calculated in Appendix C 2 

Now we can write the high temperature asymptotics for free energy 


F{T) = -e 


,T 


0.10350 + — In 


Ta 


+ 


15 c^It? 


64 2hc 65536 


+ 0{T-4 


and for entropy 


S{T) = 


e 


0.10350 + — 4 +ln^ 
64 V 2hc 


15 

65536 T‘^a? 


+ 0 {T -4 . 


(7.64) 


(7.65) 


(7.66) 


(7.67) 


Putting in these equations = 1 we arrive at the double scattering approximation for a 
perfectly conducting cylindrical shell (see eq. (7.60)). A slight distinction between the linear 
in T terms in eq. (7.60) and eq. (7.66) is due to a finite error inherent in the numerical 
methods employed in both the approaches. 

In the case of a pure dielectric cylinder (/ii = = 1, 7^ £ 2 ) the first four heat kernel 

coefficients are different from zero even in the dilute approximation [49] (small difference 
between the velocities of light inside and outside the cylinder) 


Bn = 


Qua? 12 7ra^ 2 

-4— (ci — C2) -I-5— (ci — C2) , 


Bi/2 — 


271^/4 


(ci - C2)^ 


Bi = 4- (ci - C 2 ) 


14 TT 


(ci - C2)^ 


-^ 3/2 - „2 ('21 “ C2 ) , 


16a C 2 


B2 = 0, 


Bi 


5/2 


(47r)3/2 


857 (C 1 -C 2 )' 
61440 a3 


(7.68) 


It should be noted that the coefficient B 2 vanishes only in the (ci — C 2 )^-approximation. As 
a matter of fact B 2 contains nonvanishing (ci — C 2 )^-terms and those of higher order [49]. 
Therefore the zeta regularization provides a hnite answer for the vacuum energy of a pure 
dielectric cylinder only in the (ci — C 2 )^-approximation even at zero temperature. 

The contribution to the asymptotic expansions of the first three heat kernel coefficients 
should be involved into the relevant phenomenological parameters in the general expression 
of the classical energy of a dielectric cylinder (in the same way as it has been done for a pure 
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dielectric ball). By making use of the coefficients B 3/2 and B ^/2 we get the high temperature 
asymptotics of the internal energy in the problem at hand 


U{T) = An^ 


3 T 


857 \ 

17280 T2i?2 j 


+ C>(T-2) . 


(7.69) 


where An = ni — n 2 — (c 2 — c\)/c. 

In view of considerable technical difficulties we shall not calculate the zeta function de¬ 
terminant for a pure dielectric cylinder. We recover the respective asymptotics of free energy 
by integrating the thermodynamic relation (7.13) and of entropy by using the relation (7.14). 
Pursuing this way we introduce a new constant of integration a that remains undetermined 
in our consideration 


F(T) 


-An^ 


3 T 


( , 


857 c^h^\ 
34560 T2a2 ) 


+ (P(T-2), 


(7.70) 


S{T) = An^ 


3 


1 

1 -|- a -|- In — 
he 


857 c^h^\ 
34560 T2a2 J 


+ 0 {T-^). 


(7.71) 


E. Summary of sect. VII 

We have demonstrated efficiency and universality of the high temperature expansions in 
terms of the heat kernel coefficients for the Casimir problems with spherical and cylindrical 
symmetries. All the known results in this held are reproduced in a uniform approach and in 
addition a few new asymptotics are derived (for a compact ball with ci = C 2 and for a pure 
dielectric inhnite cylinder). 

As the next step in the development of this approach one can try to retain the terms 
exponentially decreasing when T —cx). These corrections are well known, for example, for 
thermodynamic functions of electromagnetic held in the presence of perfectly conducting 
parallel plates [67] (see also eq. (7.44)). In order to reveal such terms, hrst of all the 
exponentially decreasing corrections should be retained in the asymptotic expansion (3.25) 
for the heat kernel. 

It is worth noting that in the framework of the method employed the high temperature 
asymptotics can also be constructed in the problems when the zeta regularization does not 
provide a hnite value of the vacuum energy at zero temperature, i.e., when the heat kernel 
coefficient B 2 does not vanish. Such problems were considered, for example, in ref. [190] by 
making use of the asymptotic energy densities (see also paper [191]). 

In ref. [180] it was argued that in the high temperature limit the behavior of the Casimir 
thermodynamic quantities should be the following. In the case of disjoint boundary pieces 
the free energy tends to minus infinity, the entropy approaches a constant, and the internal 
energy vanishes. Contributions to the Casimir thermodynamic quantities due to each indi¬ 
vidual connected component of the boundary exhibits logarithmic deviations in temperature 
from the behavior just described. In our consideration we were obviously dealing with an 
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individual connected component of the boundary (a sphere or cylinder). Our results corrob¬ 
orate the relevant conclusions of ref. [180] concerning the free energy and entropy. However 
the internal energy in our calculations tends to inhnity like T instead to vanish, this increase 
being caused by the respective logarithmic terms in the high temperature asymptotics of 
free energy. 

It is worth developing the methods that enable one to reveal the low temperature asymp¬ 
totics of Casimir energy by making use of the spectral zeta functions and heat kernel tech¬ 
nique. 


VIII. CONCLUSIONS 

The reality of Casimir forces and consequently the reality of relevant vacuum energy of 
quantized helds is now well founded experimentally. Therefore the task of the theory in 
this area is to develop mathematically consistent methods for calculating these quantities. 
Unfortunately we have to assert that till now one didn’t succeed in casting the dehnition of 
the Casimir energy (3.2) on a rigorous mathematical footing for boundaries of an arbitrary 
geometry. It is this reason that causes the controversies persisting here for a long time (see, 
for example, papers [61,192-197] and references therein). 

In accordance with the general concept of the quantum held theory, the subtraction 
procedure (3.2) could be rigorously specihed in the framework of proving the renormaliz- 
ability of the quantum held model under study. However it is very unlikely to implement 
this by a straightforward generalization of the renormalizability of quantum electrodynamics 
treated in unbounded Minkowski space. Point is that the boundaries or nonhomogeneous 
characteristics of the conhguration space drastically complicate the propagators [198] which 
are the kernels of the Schwinger-Dayson integral equations for complete Green’s functions. 
The proof of the renormalizability is substantially based on the analysis of these equa¬ 
tions [92,199]. 

The following consideration is also important here. The proof of renormalization theo¬ 
rem is carried out, as a matter of fact, in the framework of perturbation theory because the 
complete kernels of the Schwinger-Dayson equations are known only in the form of pertur¬ 
bation series. In order to anticipate the convergence of resulting series one has to assume 
the interaction in the theory under consideration to be small. The boundaries can also be 
treated as an interaction with appropriate classical helds. However the latter cannot be for 
sure considered as a weak one [85]. In this situation application to the calculation of the 
vacuum energy of the renormalization group technique is interesting [200,201] 

To our opinion, the present status of the Casimir calculations can be appropriately 
described by the following words due to Milton [193]: ’’Obviously we are still at the early 
stages of understanding quantum held theory. The nature of divergences in vacuum energy 
calculations is still not understood. However, there are a few established peaks that rise 
above the murky clouds of ignorance, and we should not abandon them lightly because the 
rest is obscure.” 

In this situation the spectral zeta function method and heat kernel technique are dis¬ 
tinguished here because they relay on sound mathematics, reproduce the results obtained 
in other approaches and are applicable, at least in principle, to a broad range of physical 
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problems. 

An urgent task in the Casimir studies is the development of methods that enable one 
to do calculations for boundaries without high symmetry. It is worth noting the proposal 
to employ in this held the world line formalism in QFT [202] and optical approach to the 
Casimir effect [203]. Certainly, one can anticipate here only approximate solutions. In order 
to preserve the rigorous treatment of the divergencies we suppose that it is worthy to use 
here the chain of relations between the Casimir energy, the spectral zeta function (3.13) and 
the relevant heat kernel (3.24). For the latter one can construct the equations [204] 

which provide us at least with perturbation series. The procedure of analytic continuation 
should be done here for each term of these series separately. 
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APPENDIX A: ANALYSIS OF THE DIVERGENCIES GENERATED BY lUf 

Here we reveal an important relation between linear 

and quadratic in An terms in dehned in eq. 

(5.12). 

Let us put 


I 1 

xi=y{l + —j 


I 1 

X2 = y(l-—J 


The Taylor expansion yields 

lVi(xi,X2) = Si(xi)ei(x2) - Si(xi)ei(x2) 

= -1 + (2s;e; - ste'/ - s'lei) 


Ax 


1 (s,er - .r'e,) + I w;e', - .[e'!) 


Ax = Any. 


Ax"^ 


+ 0 {Ax^ 


(Al) 


(A2) 


For brevity we have dropped the argument y of the function si and e^, and have used the 
value of the Wronskian 


W{si{y),ei{y)} = sie'i - s'lei = - 1 . 
By making use of the equation for the Riccati-Bessel functions 


(A3) 
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Wi{y) - L{l,y)wi{y) =Q , L{l,y) = l + 


/(/ + !) 
q,2 


(A4) 


we obtain 


s'l'ei - sie'” = L{l,y ), 
s'la'i-s[e'l = -L{l,y). 

Substitution of (A5) into (A2) gives 

Wi{xi.X 2 ) = -1 + [s[e'i - L{1, y)siei]Ax - ^ L{1, y) Ax^ + 0{Ax^). 
Squaring eq. (A6) one gets 


W^{xi,X 2 ) = 1 + Ai An + Bi An^ + O(An^), 


where 


A/ = y{s'lei + sie'l - 2s[e[) = 2y 
Bi = y'^L{l, y) + -Af . 


2L{l,y)siei - - (s^e;)' 


In terms of these notations we can write 


In 


AA 


= AiAn + B, 


AW? 


AA — Pi + O(An^). 


(AS) 

(A6) 

(A7) 

(AS) 

(A9) 

(AlO) 


The terms quadratic in An in eq. (AlO) exactly reproduce the function Fi{y) in eq. (9) of 
the paper [40]. It is this function that affords the whole hnite value of the Casimir energy 
in the problem under consideration. Unlike the papers [40,41,43] we didn’t introduce the 
contact terms in the dehnition of the Casimir energy and nevertheless we have reproduced 
the key function Fi{y)- It implies that the contact terms do not give a contribution into the 
hnite part of the Casimir energy in this problem. They merely cancel the terms AiAn in eq. 
(AlO). 

Now we show, without invoking the contact terms, that the Ai terms in eq. (AlO) do not 
contribute into the vacuum energy. 

Using eq. (5.18) with 6^ = 0 we introduce the notation 

OO 

W(21 + l)si{yr)ei{yp) + 1 = .-vP-P = T>{r,p,y) . (All) 

tr 


Taking into account the explicit form of the coefficients Ai dehned in eq. (AS) one can write 


An ^(2/ + 1)A; = yAn 


1=1 




-2 




+ 




dr'i drdp dp^ 


'B{r,p,y) 


+ 1 


(A12) 


r=p=l 
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When r = p = 1 the derivatives of the function T> in eq. (A12) tend to inhnity. Therefore 
a preliminary regularization should be introduced here in order to put our consideration on 
a rigorous mathematical footing. To this end we dehne the right-hand side of eq. (A12) in 
the following way 


An ^^(2/ -|- l)Ai = An lirn {V^r — 2Vrp + T^pp) 


1=1 


r = l+e/2 
p=l-e/2 


+ 1 ) 


(A13) 


where the positive constant e is a regularization parameter. From the explicit form of the 
function V{r,p,y) (see eq. (All)) it follows immediately 


lim [Vrr - Vpp) 


r=l + e/2 
p=l-e/2 


= 0 . 


(A14) 


The analogous limit for the differences 


— 'Dj.p and '^pp — '^rp 


(A15) 


also vanishes. Hence in the regularization introduced above the sum under consideration 
has the following value 


An^(2/ + l)Hz = 1. (A16) 

1=1 


It implies immediately that the term linear in An, which encounters eq. (AlO) does not 
contribute to the vacuum energy E 2 dehned in eq. (5.11). 

Now we show that the contributions to the Casimir energy given by Bi and by 
(1/4) are the same. In other words, y‘^L{l,y) in eq. (A9) does not give any hnite 

contribution to the vacuum energy. In order to prove this, we consider the expression 



(A17) 


Instead of the cutoff regularization we shall use here the analytical regularization presenting 
(A17) in the following form 


I 



y^-^dy = lim ^ 
1=1 







(AI8) 


Here the change of integration variable y = uz is done and the photon mass p is introduced. 
Further we have 
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(A19) 


I 


lim lim [(2"^+" - l)C(s - 4) - 2“^+"] 


-—lim 

24 


lim 

Aj2^0 



0 . 



In view of all this, we are left with the following scheme for calculating the Casimir 
energy of a dielectric ball in the An^-approximation. First, the An^-contribution should 
be hnd, which is given by the sum Upon changing its sign to the opposite one, 

we obtain the contribution generated by W^, when this function is in the argument of the 
logarithm. Obviously, this result would be deduced directly if one could hnd in a closed 
form the sum Yli [33]. This assertion can be represented by a symbolic formula 

In ^ ~ ^ P? + O(An^). (A20) 


The sign ~ means here the equality subject to the regularizations described above are 
employed. 


APPENDIX B: ZETA FUNCTION DETERMINANTS FOR SPHERICALLY 

SYMMETRIC BOUNDARIES 


1. Perfectly conducting sphere 


We shall use here the integral representation for the zeta function given in eq. (4.17) with 
= 1. The analytic continuation of this expression to the region Re s < 0 is performed 
by adding and subtracting from the integrand its uniform asymptotics at large u 



t{z) 


1 


(Bl) 


As a result we obtain 


C(s) — Z{s) + Cas(-s), 


(B2) 


where 

Z{s) 


a\2s sin(7rs) 
c) 2 tt 




-2s 


Cas(5) - (J 


a\ 2s 3sin(7rs) 


1=1 
OO 


/‘°° dz d f . 2 / NT 3 1 

/o I ^^{1 + ;^2)3 J ’ 


(B3) 


l-2s +8 


C 

1 /a\2s 


nz) 


4 \c 


^ poo 

/ dzz^- 

A' Jo 

s(l + s)(2 + s) [(2'+"* - 1 )Cr( 1 + 2s) - 2’+q . 


(B4) 
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When calculating C^O) one can put in eq. (B3) s = 0 everywhere except for sin(7rs), the 
latter function being substituted simply by tts. In view of this the integral in eq. (B3) is 
evaluated easy if one takes into account the limits 


\_2T{iy + l)\ 4z/2’ 4z/2 

and the asymptotics (Bl) at large As a result we obtain 


4z/2 / 4 z/2 


Differentiation of eq. (B4) with respect to s at the point s = 0 gives 

CL(0) ^ “3 + 2 2 ' 

In order to calculate the sum over I in eq. (B6) we consider an auxiliary sum 
^ r / 2 \ 2 " 

= In l-£jU^ . S(0) = 0. S(1) = Z'(0). 


where g is a parameter. Derivative of this sum with respect to q can be rewritten in the 
form 


^ ~ ^ [/ + 1/2 “ / + (l + g)/2 ^ / + 1/2 “ / + (1 -g)/2j ‘ ^ 

The summation in eq. (B9) is accomplished by making use of the following relations [112] 


^_^ 1 _ 1 

y + k X + k J X y 


+ '4){x) 


i){x + l) = i){x) +- =- 7 - 21 n 2 , 


(BIO) 


where ijj{x) is the digamma function (the Euler 'ijj function): ip{x) = {d/dx)lnr{x). This 
gives 

S'(5)=«(2-7-2ln2)-| ^([l + l])+V([|-|]) . (Bll) 


Now we integrate the both sides of eq. (Bll) over a from 0 to 1 by making use of ‘Maple’ 

S(l) = Z'(0) = 1 - I + 1 1.12 + 6cy-l). (B12) 

/,From eqs. (B2), (B7) and (B12) it follows that 

C'(0) = i-| + b.i2 + 6Cy-l)+(-5 + bn^ + ln2+2) 

= ^ In 2 + 6 Cr(- 1) + i In - = 0.38429 + ^ In -. (B13) 

o 0 ^ C ^ C 
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2. Material ball with ci = C 2 


The complete zeta function in this problem is given in eq. (4.17) (without expanding 
in powers of the parameter ). Adding and subtracting under the integral sign in this 
equation the uniform asymptotics of the integrand at large v we get 


C(s) = ( 


a\2s sin(7rs) 


c) 


271 


E"' 

i=i 


-•2s 


dz d 
dz 


z^ 


ln[l - ^ CT; {vz)] + 




4z/2 (1 + z'^Y 


+'C^Cas(’S), 


(B14) 


where the function Cas(>s) was introduced in eq. (B4). Proceeding in the same way as in the 
previous subsection we obtain for the derivative of the function (C(s) at the point s = 0 

C'(0) = ^(a+e\L(0), (B15) 

where the function S{Y^) is defined in eq. (B8). For small values of the argument we 
deduce from eq. (B8) 

1 cA 

E ^ p - 8| + OK')- (B16) 

^ 1=1 ^ 

Therefore restricting ourselves to the first order of we arrive at the final result 

C'(0) = CJO) = «' (-| + ^'“7 + '■'2 + l) ■ (B17) 


APPENDIX C: ZETA FUNCTION DETERMINANTS FOR CYLINDRICALLY 

SYMMETRIC BOUNDARIES 


1. Perfectly conducting cylindrical shell 


The spectral zeta function for this configuration has been considered in subsect. IV C. 
This function was used there for obtaining a finite value for the relevant Casimir energy. 
For this the value of the zeta function at the point s = —1/2 was calculated. Now we are 
interested in the value of C^(0)- Therefore it is convenient to represent the initial formula 
(4.46) with = 1 in a slightly different way as compared with the subsect. IV C, namely 


,2s-l 


C(^) = 


2v^c2^F(s) F(3/2 - s) 

^2s-l 


dyy^ -^\n[l - nl{y)] 


V?c2^F(s)F(3/2 




dyy^ ^*^ln[l-/i^(n|/)], 


(Cl) 


where fin{y) was defined in eq. (4.47). The first term on the right hand side of eq. (Cl) is 
an analytic function of the complex variable s in the strip —1/2< Res<l/2. Therefore 
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there is no need in analytic continuation of this expression when calculating C^O)- As 
regard to the second term in eq. (Cl) its analytic continuation to the region Re s < 0 can 
be accomplished in a standard way. We add and subtract here the uniform asymptotics of 
the integrand when n tends to inhnity 


\n[l - ^il{ny)] ~ 


y^t\y) 


+ 0{n-^), 


Ky) 


1 

yr+F' 


(C2) 


As a result we obtain 


C(^) = 


,2s-l 


dy d 


ln[l - yl{y)] 


2^/^c^^T{s)T{3/2 - s) Jq y"^^ ^ dy 

2^-1 “ dy d 


v/^c2^r(s) r(3/2 


Jo 


y 


2s-l 


dy 


ln[l - yliny)] + 


A+6 


4n^ 


,2s-l 


32^/nc■ 


2s 


(1-2s)(3-2s)Cr(2s + 1) 


r(l /2 + s) 

r(i>) ' 


(C3) 


Keeping in mind the behavior of the gamma function at the origin r(s) ~ s ^ one can easily 
find the derivative of C(-s) at the point s = 0 


1 d 

C'(0) = — / dyy— ln[l - yl{y)] 

TTO Jq dy 

Unlike the spherically symmetric boundaries, the integration is not removed in the formula 
obtained for C^O)- Therefore the first two terms in eq. (C4) can be calculated only numeri¬ 
cally 


1 

Tia 


dy \n[l - yl{y)] 


0.53490 

a 


(C5) 


Applying the FORTRAN subroutine that approximates the Bessel functions by Chebyshev’s 
polynomials we evaluate the first 30 terms in the sum in eq. (C4) 


2 ^ roo I' 

-g’7„ "0“'’'- 


Ml(ny)] + 


4n^ 


0.00554 


(C6) 


Finally gathering together all these results we have 

0.45847 3 , a 

C'O =-+ In-. 

a 32 a 2c 


(C7) 
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2. Compact infinite cylinder with ci = C 2 


Now we turn to a compact cylinder placed into unbounded medium such that the velocity 
of light is uniform on the lateral snrface of the cylinder. Proceeding as in the case of 
a cylindrical shell we rewrite the initial eqnation (4.46) in the linear approximation with 
respect to in the form 


C(^) 


,2s-l c2 


e 


2v/ic2^r(s)r(3/2 - s) io 


dyy^ '^"-^yliy) 

dy 


v^c2^r(s)r(3/2 - s) u 


d 


dyy' 


(C8) 


where y,n{y) is defined in eq. (4.47). The analytic continnation to the region Re s < 0 
is needed only for the second term in eq. (C8). Adding and snbtracting here the nniform 
asymptotics of the integrand for large n 

-yl(ny) ^ + 0{n-^), (C9) 

we obtain 


C'(o) 


-/ dyy—fiQ{y) + -^ 

TTO Jq dy na 


n Hdyy-^ 
Jo dy 


-yliny) + 


4 +6 


y"t' 

4n^ 


-;^[37-4-31n — 
32a V a 


(CIO) 


The first two terms in eq. (CIO) can again be calculated only nnmerically 


na 



— 0.28428, 
a 


(Cll) 


2^2 


Tra 


E 

n=l 


n 


dy 


-yl{ny) + 


4n^ 


= -—0.00640. 
a 


(C12) 


The final result reads 


C'(o) = 




2 r 


a 


0.28428 - 0.00640 + —( 37 - 4- 3 In — 
32 V a 


= — (0.20699 + — In — ) . 
a \ 32 2c; 


(C13) 


70 



REFERENCES 


[1] Casimir H. B. G., Proc. K. Ned. Akad. Wet. 51 (1948) 793. 

[2] Sparnaay M. J., Physica (Utrecht) 24 (1958) 751; Physics in Making, edited by A. Sar- 
lemijn and M. J. Sparnaay (North-Holland, Amsterdam) 1989. 

[3] Tabor D., Winterton R. H. S., Proc. R. Soc. London A 312 (1969) 435; Is- 
raelachvili J. N., Tabor D., Proc. R. Soc. London A 331 (1972) 19; Van Blok- 
land P. H. G. M., Overbeek J. T. G., J. Ghem. Soc. Faraday Trans. 74 (1978) 2637. 

[4] Lamoreaux S. K., Phys. Rev. Lett. 78 (1997) 5; Gomm. Mod. Phys. D; At. Mol. Phys. 
2 (2000) 247. 

[5] Mohideen U., Roy A., Phys. Rev. Lett. 81 (1998) 4549; Roy A., Mohideen U., Phys. 
Rev. Lett. 82 (1999) 4380; Harris B. M., Ghen F., Mohideen U., Phys. Rev. A 62 (2000) 
052109; Roy A., Lin G-Y, Mohideen U., Phys. Rev. D 60 (1999) 111101; Gomm. Mod. 
Phys. D; At. Mol. Phys. 2 (2000) 263. 

[6] Ghan H. B. et al., Science 291 (2001) 1941. 

[7] Bressi G., Garugno G., Onofrio R., Ruoso G., Phys. Rev. Lett. 88 (2002) 041804. 

[8] Serry F. M. et al., J. Microelectromechanical Syst. 4 (1995) 193. 

[9] Ghan H. B. et al., Phys. Rev. Lett. 87 (2001) 1801. 

[10] Ghen F., Klimchitskaya G. L., Mohideen U., Mostepanenko V. M., Phys. Rev. Lett. 
90 (2003) 160404. 

[11] Ghen F., Klimchitskaya G. L., Mohideen U., Mostepanenko V. M., Phys. Rev. A 69 
(2004) 022117. 

[12] Bimonte G., Galloni E., Esposito G., Milano L. and Rosa L., arXive: quant-ph/0406188 
v3. 

[13] Boyer T. H., Phys. Rev. 174 (1968) 1764. 

[14] Davies B., J. Math. Phys. 13 (1972) 1324. 

[15] Balian R., Duplantier B., Ann. Phys. (N.Y.) 104 (1977) 300; ibid. 112 (1978) 165. 

[16] Milton K. A., DeRaad (Jr.) L. L., Schwinger J., Ann. Phys. (N.Y.) 115 (1978) 388. 

[17] Nesterenko V. V. , Pirozhenko 1. G., Phys. Rev. D 57 (1997) 1284. 

[18] Bowers M. E., Hagen G. R., Phys. Rev. D 59 (1999) 025007. 

[19] Milton K. A., Nesterenko A. V., Nesterenko V. V., Phys. Rev. D 59 (1999) 105009. 

[20] DeRaad (Jr.) L. L., Milton K., Ann. Phys. (N.Y.) 136 (1981) 229; DeRaad (Jr.) L. L., 
Fortschr. Phys. 33 (1985) 117. 

[21] Gosdzinsky P., Romeo A., Phys. Lett. B 441 (1998) 265. 

[22] Lifshitz E. M., Zh. Eksp. Teor. Fiz. 29 (1955) 894; [Sov. Phys. JETP 2 (1956) 73]; 
Lifshitz E. M., Pitaevskii L. P., Statistical Physics (Pergamon, Oxford) 1980, Part 2. 

[23] van Kampen N. G., Nijboer B. R. A., Schram K., Phys. Lett. A 26 (1968) 307. 

[24] Schwinger J., DeRaad (Jr.) L. L. and Milton K. A., Ann. Phys. (N.Y.) 115 (1978) 1. 

[25] Brevik L, Nesterenko V. V. and Pirozhenko 1. G., J. Phys. A 31 (1998) 8661. 

[26] Nesterenko V. V. and Pirozhenko 1. G., Phys. Rev. D 60 (1999) 125007. 

[27] Brevik L, Kolbenstvedt H., Ann. Phys. (N.Y.) 143 (1982) 179; ibid. 149 (1983) 237; 
Gan. J. Phys. 62 1984 805; ibid. 63 (1985) 1409. 

[28] Brevik 1. and Nyland G. H., Ann. Phys. (N.Y.) 230 (1994) 321. 

[29] Brevik 1. and Einevoll G., Phys. Rev. D 37 (1988) 2977. 

[30] Brevik 1. and Glausen L, Phys. Rev. D 39 (1989) 603. 


71 



[31] Brevik I. J. Phys. A 20 (1987) 5189. 

[32] Brevik I. and Yousef T. A., J. Phys. A 33 (2000) 5819. 

[33] Klich L, Phys. Rev. D 61 (2000) 025004. 

[34] Klich L, Phys. Rev. D 64 (2001) 045001. 

[35] Klich I. and Romeo A., Phys. Lett. B 476 (2000) 369. 

[36] Klich I., Feinberg J., Mann A. and Revzen M., Phys. Rev. D 62 (2000) 045017. 

[37] Milton K. A., Ann. Phys. (N.Y.) 127 (1980) 49. 

[38] Milton K. A. and Ng Y. J., Phys. Rev. E 57 (1998) 5504. 

[39] Barton G., J. Phys. A 32 (1999) 525. 

[40] Brevik I., Marachevsky V. N. and Milton K. A., Phys. Rev. Lett. 82 (1999) 3948. 

[41] Brevik 1. and Marachevsky V. N., Phys. Rev. D 60 (1999) 085006. 

[42] Lambiase G., Scarpetta G. and Nesterenko V. V., Mod. Phys. Lett. 16 (2001) 1983. 

[43] Marachevsky V. N., arXive: hep-th/9909210, Preprint SPb-IP-99-11; Phys. Scripta 64 
(2001) 2005. 

[44] Hpye J. S. and Brevik L, J. Stat. Phys. 100 (2000) 223. 

[45] Hpye J. S., Brevik 1. and Aarseth J. B. Phys. Rev. E 63 (2001) 051101. 

[46] Nesterenko V. V., Lambiase G. and Scarpetta G., Phys. Rev. D 64 (2001) 025013. 

[47] Barton G., Phys. Rev. A 64 (2001) 032103. 

[48] Barton G., J. Phys., A 34 (2001) 4083. 

[49] Bordag M. and Pirozhenko 1. G., Phys. Rev., D 64 (2001) 025019. 

[50] Gavero-Pelaez 1. and Milton K. A., arXiv: hep-th/0412135. 

[51] Romeo A. and Milton K. A., arXive: hep-th/0504207. 

[52] Derjaguin B. V. and Abriksova 1. L, Sov. Phys. JETP 3 (1957) 819; Derjaguin B. V., 
Sci. Am. 203 (1960) 47. 

[53] Gies H., Langfeld K. and Moyaerts L., JHEP 0306 (2003) 018. 

[54] Schaden M. and Spruch L., 

[55] Saharian A. A., Phys. Rev. D 63 (2001) 125007; Int. J. Mod. Phys. A 19 (2004) 4301. 

[56] Saharian A. A., IGTP preprint, IG/2000/14; hep-th/0002239. 

[57] Ahmedov H. and Burn 1. H., J. Math. Phys. 44 (2003) 5487. 

[58] Jaffe R. L. and Scardicchio A.,e arXiv: hep-th/0501171. 

[59] Milton K. A., J. Phys. A 37 (2004) R209. Phys. Rev. Lett. 84 (2000) 459. 

[60] Bordag M., Mohideen U. and Mostepanenko V. M. Phys. Rep. 353 (2001) 1. 

[61] Milton K. A., The Casimir Effect: Physical Manifestations of Zero-Point Energy, 
(World Scientific, Singapoore) 2001. 

[62] Mostepanenko V. M. and Trunov N. N., The Casimir effect and its applications 
(Glarendon, Oxford) 1997. 

[63] Grib A. A., Mamaev S. G. and Mostepanenko V. M., Vacuum Quantum Effects in 
Strong Eields (Friedman Laboratory Publishing, St. Petersburg) 1994. 

[64] Belifante F. J., Am. J. Phys. 55 (1987)134; Elizalde E. and Romeo A., Am. J. Phys. 
59 (1991) 711. 

[65] Milonni P. W., The Quantum Vacuum: An Introduction to Quantum Electrodynamics 
(Academic Press, San Diego) 1994. 

[66] Kardar M. and Golestanian R., Rev. Mod. Phys. 71 (1999) 1233. 

[67] Plunien G., Muller B. and Greiner W., Phys. Rep. 134 (1986) 87. 

[68] Roberts M. D., Vacuum Energy, arXive: hep-th/0012062. 


72 



[69] Leipzig 1995, Quantum Field Theory Under the Influence of External Conditions, Pro¬ 
ceedings of the 3d Workshop on Quantum Field Theory under the Influence of External 
Conditions (Leipzig, Germany, 18 - 22 Sep 1995), edited by M. Bordag (Teubner, Stut- 
gart) 1996. 

[70] The Casimir effect 50 years later, Proceedings of the 4th Workshop on Quantum Field 
Theory under the Influence of External Conditions (Leipzig, Germany, 14 - 18 Sep 
1998), edited by M. Bordag (World Scientihc, Singapore) 1999. 

[71] Proceedings of the 5th Workshop on Quantum Field Theory under the Influence of 
External Conditions (Leipzig, Germany, 10-14 Sep 2001), edited by M. Bordag Int. J. 
Mod. Phys. A 17 (2002) No. 6 and 7. 

[72] Norman 2003, Quantum Field Theory Under the Influence of External Conditions, 
Proceedings of the 6th Workshop on Quantum Field Theory under the Influence of 
External Conditions (Norman, Oklahoma, USA, 15-19 Sep 2003), edited by K. A. Mil- 
ton (Rinton Press, Paramus, NJ, USA) 2004. 

[73] Planck M., Ann. d. Phys. 37 (1912) 642; Milonni P. W. and Shih M.-L., Amer. J. 
Phys. 59 (1991) 684. 

[74] Einstein A., Stern O., Ann. d. Phys. 40 (1913) 551. 

[75] Debye P., Ann. d. Phys. 43 (1914) 49. 

[76] Mulliken R. S., Nature 114 (1924) 349. 

[77] Einstein A., Ann. d. Phys. 22 (1907) 180. 

[78] Debye P., Ann. d. Phys. 39 (1912) 789. 

[79] Matthews P. T. and Salam A., Phys. Rev. 94 (1954) 185. 

[80] Gunn J. C., Rep. Prog. Phys. 18 (1955) 127. 

[81] Dirac P. A. M., Lectures on Quantum Field Theory (Published by Belfer Graduate 
School of Science, Yeshiva University, New York) 1967. 

[82] Schwinger J., Lett. Math. Phys. 1 (1975) 43; ibid. 24 (1992) 59; ibid. 24 (1992) 227. 

[83] Milonni P. W., Phys. Rev. A 25 (1982) 1315. 

[84] Milonni P. W., Phys. Scripta T 21 (1988) 102. 

[85] Jaffe R. L. preprint MIT-CTP-3614; arXiv: hep-th/0503158. 

[86] Barash Yu. S., Van der Waals Forces (Nauka, Moscow) 1988. 

[87] Hertz Heinrich, Die Principien der Mechanik in neuem Zusammenhange dargestellt, 
Gesam. Werke, Bd. HI (Leipzig) 1910. 

[88] Sommerfeld A., Vorlesungen iiber Theoretische PhysikBd. I, Mechanik, vierte Auflage 
(Dieterich’sche Verlagsbuchhandlung, Wiesbaden) 1949. 

[89] Marshall T. W., Proc. R. Soc. Lond. A 276 (1963) 475. 

[90] Boyer T. H., Phys. Rev. 182 (1969) 1374. 

[91] Marshall T. W., Pros. Camb. Philos. Soc. 61 (1965) 537; Nuovo Cimento 38 (1965) 
206; Henry L. L. and Marshall T. W., Nuovo Cimento 41 (1966) 188; Boyer T. H., 
Phys. Rev. 186 (1969) 1304; Phys. Rev. D 1 (1970) 1526; ibid. 1 (1970) 2257. 

[92] Itzykson C. and Zuber J.-B., Quantum Field Theory (McGraw-Hill, New York) 1980. 

[93] Birrell N. D. and Davies P. C. W., Quantum Fields in Curved Space (Cambridge 
University Press, Cambridge, UK) 1982. 

[94] Kirsten K., Spectral Functions in Mathematics and Physics (Chapman & Hall/CRC, 
Boca Raton, FL, USA) 2001. 

[95] Elizalde E., Odintsov S. D., Romeo A., Bytsenko A. A. and Zerbini S.,, Zeta Regular- 


Id 



ization Technique with Applications (World Scientific, Singapore) 1994. 

[96] Elizalde E., Ten Physical Applications of Spectral Zeta Functions (Springer, Berlin) 
1995. 

[97] Vassilevich D. V., Pliys. Kept. 388 (2003) 279. 

[98] Santangelo E. M., Theor. Math. Phys. 131 (2002) 527; hep-th/0104025. 

[99] Blau S. K., Visser M. and Wipf A., Nucl. Phys. B 310 (1988) 163. 

[100] Bytsenko A. A., Cognola G., Vanzo L. and Zerbini S. Phys. Rep. 266 (1996) 1. 

[101] Bytsenko A. A., Cognola G., Elizalde E., Moretti V. and Zerbini S. Analytic Aspects 
of Quantum Fields (World Scientihc, Singapore) 2003. 

[102] Elizalde E., Odintsov S. D. and Romeo A. J. Math. Phys. 37 (1996) 1128. 

[103] Dowker J. S. and Critchley R. Phys. Rev. D 13 (1976) 3224. 

[104] Hawking S. W., Commun. Math. Phys. 55 (1977) 133. 

[105] Whittaker E. T. and Watson G. N., A Course of Modern Analysis (Cambridge Uni¬ 
versity Press) 1969. 

[106] Cognola G, Vanzo L. and Zerbini S., Phys. Rev. D 52 (1995) 4548; Zerbini S., Con- 
gola G. and Vanzo L., Phys. Rev. D 54 (1996) 2699; Cheeger J., J. Diff. Geom. 18 
(1983) 576. 

[107] Milton k., Odintsov S. D. and Zerbini S., Phys. Rev., D 65 (2002) 065012. 

[108] Frolov V. P. and Serebriany E. M., Phys. Rev., D 35 (1987) 3779. 

[109] Nesterenko V. V., Lambiase G. and Scarpetta G., Ann. Phys. (N.Y.) 298 (2002) 403. 

[110] Cognola G., Vanzo L. and Zerbini S., J. Math. Phys.33 (1992) 222. 

[111] Beneventano C. G. and Santangelo E. M., Int. J. Mod. Phys. A 11 (1996) 2871. 

[112] Gradshteyn I. S. and Ryzhik I. M., Table of Integrals, Series, and Products, 6th ed. 
(Academic Press, New York) 2000. 

[113] Moss I. G., Class. Quantum Grav. 6 (1989) 759; Moss I. G. and Poletti S. J., Phys. 
Lett. B 333 (1994) 326. 

[114] D’Eath P. D. and Esposito G. V. M., Phys. Rev. D 43 (1991) 3234; D’Eath P. D. 
and Esposito G. V. M., Phys. Rev. D 44 (1991) 1713; Esposito G., Quantum Gravity, 
Quantum Cosmology and Lorentzian Geometries (Springer, Berlin) 1994. 

[115] Esposito G. and Kirsten K., Phys. Rev., D 66 (2002) 085014. 

[116] Nesterenko V. V., Pirozhenko I. G. and Dittrich J., Class. Quantum Grav. 20 (2003) 
431. 

[117] See the article due to Milton K. A. in ref. [70] pp. 20 - 36. 

[118] Nesterenko V. V. and Pirozhenko I. G., Phys. Rev. D 57 (1998) 1284. 

[119] Borgnis E., Papas C. H., Electromagnetic Waveguides and Resonators, In Encyclopedia 
of Physics, edited by S. Fliige (Springer Verlag, Berlin) 1958; v. XVI, pp. 285-422. 

[120] Jackson J. D., Classical Electrodynamics, 2nd edn (Wiley, New York) 1975. 

[121] Rubinstein I., Rubinstein L. Partial Differential Equations in Classical Mathematical 
Physics (Cambridge University Press, Cambridge, UK) 1998. 

[122] Stratton J. A., Electromagnetic Theory, (McGraw-Hill, New York) 1941. 

[123] Debye P., Ann. d. Phys. 30 (1909) 57. 

[124] Lambiase G., Nesterenko V. V. and Bordag M., J. Math. Phys. 40 (1999) 6254. 

[125] Adler S. and Piran T., Rev. Mod. Phys. 56 (1984) 1. 

[126] Fishbane P., Gasiorowics S. and Kaus P., Phys. Rev. D 36 (1987) 251; ibid. 37 (1988) 
2623. 


74 



[127] Chodos A., Jaffe R. L., Thorn C. B. and Weisskopf V., Phys. Rev. D 9 (1974) 3471; 
Hasenfratz P. and Kuti J., Phys. Rep. C 40 (1978) 75; Hiller J. R., Ann. Phys. (N.Y.) 
144 (1982) 58; Phys. Rev. D 30 (1984) 1520. 

[128] Milton K. A., Phys. Rev. D 22 (1980) 1441; ibid. 22 (1980) 1444. 

[129] Bender C. M. and Hays P., Phys. Rev. D 14 (1976) 2622. 

[130] Bordag M., Elizalde E., Kirsten K. and Leseduarte S., Phys. Rev. D 56 (1997) 4896. 

[131] Elizalde E., Bordag M. and Kirsten K., J. Phys. A 31 (1998) 1743. 

[132] Landau L. D. and Peierls R., Zs. Physik 62 (1930) 188. 

[133] Bordag M., J. Phys. A 28 (1995) 755. 

[134] Elizalde E., Leseduarte S. and Romeo A., J. Phys. A 26 (1993) 2409. 

[135] Leseduarte S. and Romeo A., J. Phys. A 27 (1994) 2483. 

[136] Magnus W., Oberhettinger F. and Soni R. P., Formulas and Theorems for the Special 
Functions of Mathematical Physics, 3rd ed. (Springer Verlag, New York) 1966. 

[137] Olver F. W. J., Phil. Trans. Royal Soc. London A 247 (1954) 328. 

[138] Abramowitz M. and Stegun 1. A. (eds) Handbook of Mathematical Functions (New 
York, Dover) 1972. 

[139] Casimir H. B. G., Physica 19 (1956) 846. 

[140] Milton K. A. and Ng Y. J., Phys. Rev. E 55 (1997) 4207. 

[141] Nesterenko V. V. and Pirozhenko 1. G., J. Math. Phys. 41 (2000) 4521. 

[142] Leseduarte S. and Romeo A., Ann. Phys. (N.Y.) 250 (1996) 448. 

[143] Sen S., J. Math. Phys. 22 (1981) 2868; ibid. 25 (1984) 2000. 

[144] Bender G. M. and Milton K. A., Phys. Rev. D 50 (1994) 6547; Milton K. A., Phys. 
Rev. D 55 (1997) 4940. 

[145] Gognola G., Elizalde E. and Kirsten K., J. Phys. A 34 (2001) 7311. 

[146] Bordag M., Kirsten K. and Vassilevich D., Phys. Rev. D 59 (1999) 085011. 

[147] Milonni P. W., Schaden M. and Spruch L., Phys. Rev. A 59 (1999) 4259. 

[148] Barton G., J. Phys. A 34 (2001) 4083; ibid. 34 (2001) 5781. 

[149] Schwinger J., Proc. Natl. Acad. Sci. USA, 89 (1992) 4091, 11118; 90 (1993) 958, 2105, 
4505, 7285; 91 (1994) 6473. 

[150] Bordag M., Nesterenko V. V. and Pirozhenko 1. G., Phys. Rev. D 65 (2002) 045011. 

[151] Phillips M., Classical Electrodynamics, in Encyclopedia of Physics, edited by S. Fliige 
(Springer Verlag, Berlin) 1962, v. IV, p. 1-108. 

[152] Heyn E., Math. Nachrichten 13 (1955) 25. 

[153] Nesterenko V. V., Lambiase G. and Scarpetta G., J. Math. Phys. 42 (2001) 1974. 

[154] Dowker J. S., preprint MUTP/2000/2; arXive: hep-th/0006138. 

[155] Gandelas R, Ann. Phys. (N.Y.) 143 (1982) 241; ibid. 167 (1986) 257; 

[156] Ford L. H. and Svaiter N. F., Phys. Rev. D 58 (1998) 065007. 

[157] Sommerfeld A., Proc. London Math. Soc. 28 (1897) 417; Garslaw H. S. Proc. London 
Math. Soc. 20 (1898) 121; Fedosov B. V., Sov. Mat. Dokl. 4 (1963) 1092; ibid. 5 
(1964) 988; Dowker J. S., J. Phys. A 10 (1977) 115; Phys. Rev. D 18 (1978) 1856; 
Gheeger J., J. Diff. Geom. 18 (1983) 575. 

[158] Howls G. J.and Trasler S. A., J. Phys. A 31 (1998) 1911. 

[159] Gognola G., Kirsten K. and Vanzo L., Phys. Rev. D 49 (1994) 1029; Zerbini S., 
Gognola G. and Vanzo L., Phys. Rev. D 54 (1996) 2699; Fursaev D. V., Phys. Lett. B 
334 (1994) 53; Glass. Quantum Grav 11 (1994) 1431; ibid. 14 (1997) 1059. 


75 



[160] Frolov V. P. and Serebriany E. M., Phys. Rev. D 35 (1987) 3779. 

[161] Neterenko V. V., Lambiase G. and Scarpetta G., Ann. Phys. (N.Y.) 298 (2002) 403. 

[162] Waechter R. T., Proc. Gamb. Phil. Soc. 72 (1972) 439. 

[163] Dowker J. S., Glass. Quantum Grav. 11 (1994) 557; Gommun. Math. Phys. 162 (1994) 
633. 

[164] Lukosz W., Physica 56 (1971) 109; Z. Phys. 258 (1973) 99; ihid. 262 (1973) 327. 

[165] Ruggiero J. R., Zimerman A. H. and Villani A., Rev. Bras. Fis. 7 (1977) 663; J. Phys. A 
13 (1980) 761. 

[166] Mamaev S. G. and Trunov V. V., Theor. Math. Phys. 38 (1979) 228. 

[167] Ambjprn J. and Wolfram S., Ann. Phys. (N.Y.) 47 (1983) 1. 

[168] Garuso F., Neto N. P., Svaiter B. F. and Svaiter N. F., Phys. Rev. D 43 (1991) 1300. 

[169] Garuso F., De Paola R. and Svaiter N. F., Int. J. Mod. Phys. A 14 (1999) 2077. 

[170] Actor A. A., Ann. Phys. (N.Y.) 230 (1994) 303. 

[171] Actor A. A. and Bender I., Phys. Rev. D 52 (1995) 3581. 

[172] Li X., Gheng H. and Zhai X., Phys. Rev. D 56 (1997) 2155. 

[173] Queiroz H., da Silva J. G., Khanna F. G., Revzen M. and Santana A. E., arXive: 
hep-th/0311246. 

[174] Ahmedov H. and Burn I. H., J. Math. Phys. 45 (2004) 965. 

[175] Ahmedov H. and Burn I. H., J. Math. Phys. 46 (2005) 022303; ibid. 022304. 

[176] Mehra J., Physica 37 (1967) 145. 

[177] Brown L.S. and Maclay G. J., Phys. Rev. 184 (1969) 1272. 

[178] Francia M. Be., Phys. Rev. B, 50 (1994) 2908. 

[179] Brevik I., Milton K., Nojiri S. and Odintsov S. B., Nucl. Phys., B599 (2001) 305; 
Brevik I., Milton K., and Odintsov S. B., Ann. Phys. (N.Y.), 302 (2002) 120. 

[180] Feinberg J., Mann M. and Revzen M., Ann. Phys. (N.Y.), 288 (2001) 103. 

[181] Kapusta J. L., Finite temperature field theory (Gambridge University Press, Gambridge, 
England) 1992. 

[182] Feynman R. P., Statistical Mechanics (W. A. Benjamin, Inc. Advanced Book Program, 
Reading, Massachusetts) 1972. 

[183] Bowker J. S. and Kennedy G., J. Phys. A 11 (1978) 895. 

[184] Bordag M., Kirsten K. and Elizalde E., J. Math. Phys. 37 (1996) 895. 

[185] Particle Bata Group, Groom B. et ah, Eur. Phys. J. C, 15, No. 1&4 (2000) 73. 

[186] Bordag M., Nesterenko V. V. and Pirozhenko I. G., Nucl. Phys. Proc. Suppl. 104 
(2002) 228. 

[187] Bordag M., Geyer B., Kirsten K. and Elizalde E., Gommun. Math. Phys. 179 (1996) 
215. 

[188] Bordag M., Kirsten K. and Bowker S., Gommun. Math. Phys. 182 (1996) 371. 

[189] See the contribution by Barton G. in ref. [71]. Marachevsky V. N., Int. J. Mod. Phys. 
A 17 (2002) 786; Rezaeian A. H., Saharian A. A., Glass. Quantum Grav. 19 (2002) 
3625. 

[190] Balian R. and Bloch G., Ann. Phys. (N.Y.) 60 (1970) 401; ibid. 63 (1971) 592; ibid. 
64 (1971) 271; ibid. 69 (1972) 76. 

[191] Be Francia M., Phys. Rev. B 50 (1994) 2908. 

[192] Milton K. A. Phys. Rev. B 68 (2003) 065020. 

[193] Milton K. A., Finite Casimir Energies in Renormalizable Quantum Field Theory, In- 


76 



vited talk given at ‘Marcel Grossmann X’ (Rio de Janeiro, July 2003); hep-th/0401117. 

[194] Graham N., Jaffe R. L., Khemani V., Quandt M., Scandurra M. and Weigel H., Phys. 
Lett. B 572 (2003) 196. 

[195] Brevik L, Jensen B. and Milton K. A., Phys. Rev. D 64 (2001) 088701. 

[196] Hagen G. R., Phys. Rev. D 61 (2000) 065005; Eur. Phys. J. G 19 (2001) 677; arXive: 
hep-th/0004079. 

[197] Fulling S. A., J. Phys. A 36 (2003) 6857. 

[198] Bordag M., Robaschik D. and Wieczorek E., Ann. Phys. (N.Y.) 165 (1985) 192. 

[199] Bogoliubov N. N., Shirkov D. V., Introduction to the Quantum Fields (Interscience, 
New York) 1959. 

[200] Brevik L, Granda L. and Odintsov S. D., Phys. Lett., B 367 (1996) 206. 

[201] Gherednikov 1. O., Acta. Phys. Slov., 52 (2002) 221; arXive: hep-th/0206245. 

[202] Gies H., Langfeld K., Moyaerts L., JHEP 0306 (2003) 018. 

[203] Jaffe R. L., Scardicchio A., Phys. Rev. Lett. 92 (2004) 070402. 

[204] Pirozhenko 1. G., Nesterenko V. V., Bordag M., J. Math. Phys., 46 (2005) 042305. 


77 



